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^ ABSTRACT 

In this second paper in a series we present measurements of spectral features of 432 
i 1 low-redshift (z < 0.1) optical spectra of 261 Type la supernovae (SNe la) within 20 d 

of maximum brightness. The data were obtained from 1989 through the end of 2008 as 
r \ part of the Berkeley SN la Program (BSNIP) and are presented in BSNIP I (Silverman 

et al., submitted). We describe in detail our method of automated, robust spectral fea- 
i-C ture definition and measurement which expands upon similar previous studies. Using 

this procedure, we attempt to measure expansion velocities, pseudo-equivalent widths 
q (pEW), spectral feature depths, and fluxes at the centre and endpoints of each of 

Jh nine major spectral feature complexes. A sanity check of the consistency of our mea- 

"o^ surements is performed using our data (as well as a separate spectral dataset). We 

investigate how velocity and pEW evolve with time and how they correlate with each 

other. Various spectral classification schemes are employed and quantitative spectral 
r __l differences among the subclasses are investigated. Several ratios of pEW values are 

calculated and studied. The so-called Si II ratio, often used as a luminosity indicator 
Qs, (Nugent et al. 1995), is found to be well correlated with the so-called "SiFe" ratio 

{NJ and anticorrelated with the analogous "SSi ratio." Furthermore, SNe la that show 

t-H strong evidence for interaction with circumstellar material or an aspherical explosion 

04 are found to have the largest near-maximum expansion velocities and pEWs, possibly 

linking extreme values of spectral observables with specific progenitor or explosion sce- 

narios. A discussion of the relative merits of various classification schemes is presented, 
{N) and we find that purely spectroscopic classification schemes are useful in identifying 

i— I the most peculiar SNe la. However, in almost all spectral parameters investigated 

the full sample of objects spans a nearly continuous range of values. Comparisons to 
. £h previously published theoretical models of SNe la are made and some of the predic- 

tions of these models match the observations presented here. We conclude with a brief 
S— i discussion of how the measurements performed herein and the possible correlations 

presented will be important for future SN surveys. 
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- cosmology:observations - distance scale 



1 INTRODUCTION 

Type la supernovae (SNe la) have been particularly useful 
in recent years as a way to accurately measure cosmolog- 
ical parameters (e. g., |Astier et aL 2006 Riess et al. 2007} 



2009 



Wood-Vasey et al.||2007| |Hicken et al.||2009| |Kessler et al.| 



Amanullah et al.|2010| |Suzuki et al.|2011 1 , and led to 



the discovery of the accelerating expansion of the Universe 
(Riess et al. 1998 Perlmutter et al. 19991. Broadly speak- 
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ing, SNe la are the result of thermonuclear explosions of 
C/O white dwarfs (WDs) (e.g.,|Hoyle fc Fowler|[l96"o"l [Col- 



gate fc McKee| 1969| [Nomoto et al.|1984| see |Hillebrandt fc 



Niemeyer 20001 for a review). However, we still lack a de- 
tailed understanding of the progenitor systems and explo- 
sion mechanisms, as well as how differences in initial condi- 
tions create the variance in observed properties of SNe la. To 
solve these problems, and others, detailed and self-consistent 
observations of many hundreds of SNe la are required. 

The cosmological application of SNe la as precise dis- 
tance indicators relies on being able to standardise their lu- 
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minosity. Phillips ( 1993 1 showed that the light-curve decline 
rate is well correlated with luminosity at peak brightness for 
most SNe la, the so-called "Phillips relation." However, this 
simple empirical relation relies on photometry alone, and 
it may be possible to refine the relation with the addition 
of spectral observations. Many comparisons of spectral fea- 
tures and studies of the temporal evolution of these features 
in low-redshift SN la have been performed in the past (e.g., 



Barbon et al.||1990 |Branch & van den Bergh 


1993| Nugent 


et al.|1995||Hatano et al.|2000||Folatelli|2004| 


Benetti et al. 


2005[ |Bongard et al.||2006| |Hachinger et al. |2006| |Bronder 


et al.|2008 Foley et al.|2008 IBranch et al.|2009||Wang et al. 


2009 Walker et al.|2011| INordin ct al. 2011b 


Blondin et al. 



there has been similar work with SNe la at higher redshifts 



(e-g-, 


Hook et al.||2005| 


Blondin et al. 2006 Altavilla et al. 


2009 


Garavini et al. 


2007 


|Bronder et al. 2008; Walker et al. 


2011 


Konishi et al. 


2011 


. Many of these studies aimed to 



find a "second parameter" in SN la spectra which would 
make our measurements of the distances to SNe la even 
more precise. 

However, most of these studies utilised relatively small 
and heterogeneous datasetsj]] Using the self-consistently ob- 
served and reduced low-redshift (z ?C 0.2) optical SN la spec- 
tra from the Berkeley Supernova la Program (BSNIP; Sil- 
verman et al., submitted), we can accurately and robustly 
measure various spectral features. These measurements can 
then be used to investigate how the spectral observables cor- 
relate with each other and with the objects' previously de- 
termined spectral subclasses based on different classification 
schemes. 

We provide an overview of the dataset used for this anal- 
ysis in Section [2] and we describe in detail our automated 
and robust procedure for measuring multiple aspects of each 
spectral feature in Section [3] Our resulting measurements 
are described in Section [4] Section [5] presents the tempo- 
ral evolution of these measured values, and how they corre- 
late with each other and with previously determined spec- 
tral classifications. We discuss our conclusions in Section [6j 
specifically summarising the main results from our analy- 
sis in Section |6.1| Finally, we attempt to answer questions 
regarding how SNe la should be spectroscopically classified 



(Section 6.2|, whether theoretical models can explain the 



spectra of SNe la and the correlations we find (Section 6.3 1, 
and how our analysis of spectral features will be important 
for future SN surveys (Section |6.4[ ). Forthcoming BSNIP pa- 
pers will utilise the spectral measurements described here 
and examine the correlations between these and other ob- 
servables (such as photometry and host-galaxy properties) . 



2 SPECTRAL DATASET 

The SN la spectra that are used in this study all come from 
BSNIP and are published in BSNIP I (Silverman et al., sub- 
mitted) . The majority of the spectra were obtained with the 
Shane 3 m telescope at Lick Observatory using the Kast dou- 
ble spectrograph ( |Miller fc Stone||1993] ). The typical wave- 
length coverage is 3300-10,400 A with resolutions of ~11 



1 The significant exception to this is the study by |Blondin et ah] 
( 2011ft. 



and ~6 A on the red and blue sides, respectively (crossover 
wavelength ~ 5500 A) . The relative spectrophotometry of 
the BSNIP data is quite accurate, though the absolute spec- 
trophometry is only correct in a handful of spectra. Thus, 
flux measurements alone may be inaccurate, but ratios of 
flux values should be quite precise. Some of the spectra 
examined here have had residual host-galaxy contamina- 
tion removed using our "colour matching" technique. This 
method uses photometry of the host galaxy of a SN to cor- 
rect for any contamination that remains after our normal 
data-reduction procedure (which often removes the major- 
ity of host-galaxy light). Further information regarding the 
observations, data reduction, spectrophotometric accuracy, 
and host-galaxy corrections can be found in BSNIP I. 

For this study we required that a spectrum be within 
20 d (rest frame) of maximum brightness, using the red- 
shift and Julian Date of maximum presented in Table 1 
of BSNIP I. The only SNe which we ignored a priori 
were the extremely peculiar SN 2000cx (e.g 



20011, SN 2002cx (e.g., Li et al. 2003 Jha et al 



and SN 2008ha (e.g., Foley et al.|2009 Valenti et al 



Li et al. 



SN 2005hk (e.g., Chornock et al.|2006 Phillips et al 



2006), 



2007), 



2009) 



After removing these objects, we were left with 458 spectra 
(147 of which were corrected for host-galaxy contamination) 
of 271 SNe la, and we attempted to measure their spectral 
features. 

Of these, there were some spectra which did not pass our 



minimum signal-to-noise ratio (S/N) cut (see Section 3.31, 
or whose wavelength range did not sufficiently cover any of 



the features we wanted to measure (see Section 3.1 1. In ad- 



dition, the endpoints of each feature sometimes fell outside 



our allowed boundaries (see Section 3.2|, or every spectral 



feature that was measured was deemed to have a poorly de- 



fined continuum (see Section 3.6 1. After removing these data 
there remain 432 spectra of 261 SNe la with a "good" fit for 
at least one spectral feature. The largest redshift of these 
observations is ~0.1 (even though the full BSNIP dataset 
contains SNe with 0.1 < z < 0.2). The earliest spectrum 
successfully fit in this work has a rest-frame age of about 
— 12.7 d. A summary of these SNe la, their ages, and spec- 
tral classifications based on various classification schemes 
can be found in Appendix [A"] 



3 MEASUREMENT PROCEDURE 
3.1 Measured Features 

Previous studies similar to this one have split optical SN la 
spectra near maximum brightness into eight or nine major 



absorption feature complexes (e.g., Riess et al. 1997 Fo- 
latelli|2004| |Hachinger et al.|2006[ ). All of these are features 
are actually blends of multiple spectral transitions, but each 
absorption complex itself is often distinct enough from the 
others for its properties to be measured independently. We 
follow previous studies' naming convention for the measured 
features by referring to each one by an ion or spectral line 
responsible for the majority of the absorption. The nine fea- 
tures we attempt to fit in each observation are labelled on a 
spectrum of the "normal" Type la SN 2002ha (taken 1 d be- 
fore maximum brightness) in Figure [I] Each feature's name 
and reference number, along with its rest wavelength (used 
to determine expansion velocities), is presented in Table [T] 
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Table 1. Spectral Features and Boundaries 



Feature Name 


Feature # 


Rest Wavelength a (A) 


Blue Boundary 1 " (A) 


Red Boundary b (A) 


Ca II H&K 


fl 


3945.28 


3400-3800 


3800-4100 


Si II A4000 


f2 


4129.73 


3850-4000 


4000-4150 


Mg II 


f3 


C 


4000-4150 


4350-4700 


Fe II 


fi 


. . , c 


4350-4700 


5050-5550 


S II "W" d 


f5 


5624.32 


5100-5300 


5450-5700 


Si II A5972 


ffi 


5971.85 


5400-5700 


5750-6000 


Si II A6355 


f7 


6355.21 


5750-6060 


6200-6600 


O I Triplet 


f8 


7773.37 


6800-7450 


7600-8000 


Ca II Near-IR Triplet 


f9 


8578.75 


7500-8100 


8200-8900 



a The rest wavelengths are weighted averages of the strongest spectral lines that give rise to each absorption feature. 
b These boundaries are necessary in order to account for variations in spectral feature width and expansion velocity 
among SNe, as well as the temporal evolution of these values. 

c This feature is a blend of so many spectral lines that a single reference wavelength is practically meaningless (and 
thus an expansion velocity cannot be accurately determined). 

d The two broad absorptions that make up the S II "W" are fit using a single spline, but we calculate the expansion 
velocity of the absorption complex using the minimum of the redder of the two features relative to its rest wavelength. 



06 
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Figure 1. The nine features we attempt to fit in each observation 
shown on a spectrum of the "normal" Type la SN 2002ha, taken 
1 d before maximum brightness, from BSNIP I. The spectrum has 
had its host-galaxy recession velocity removed and has been cor- 
rected for Milky Way reddening according to the values presented 
in Table 1 of BSNIP I and assuming that the extinction follows 
the|Carde!li et al.| (|1989} extinction law modified by |0'Donnell| 
< |1994| |. See Table|l|for more information regarding each feature. 



There have been a significant number of spectroscopi- 
cally peculiar SNe la observed whose spectra near maximum 
brightness do not contain some of the nine features or ex- 
hibit extremely weak absorptions from certain features (for 
a review of these objects, see [Filippenko| ,1997), We con- 
sider an object "spectroscopically normal" if it is classified 
as "la-norm" by the SuperNova IDentification code (SNID; 
|Blondin fe Tonry|[2007| as implemented in BSNIP I. In the 
current study we have 213 la-norm objects and 5 "la" ob- 
jects for which we were unable to determine a definitive 
subtype in BSN IP I. There are 24 SN 1991bg-like objects 
("Ia-91bg," e.g., |Filippenko et al.||1992b| |Leibundgut et aF 



1993[ ) which represent the usually underluminous SNe la. 
We also have 6 SN 1991T-like objects ( "Ia-91T," e.g., |Filip- 
penko et al^|1992a| [Phillips et al.||1992|> and 13 SN 1999aa- 
like objects ("Ia-99aa7' e. g., |Li et al.|[2001| |Strolger et"aT] 
2002 Garavini et al.|2004 1 which together represent the of- 
ten overluminous SNe la. These classifications are listed in 
the "SNID (Sub)Type" column of Table [AT] See BSNIP I 
for more information regarding our implementation of SNID 
and the various spectroscopic subtype classifications. We in- 
clude these peculiar objects in our study in an attempt to 
better quantify the extent of the spectral peculiarities among 
these types of SNe la. However, as mentioned in Section [2] 
there are some objects which are so spectroscopically pecu- 
liar (compared with "normal" SNe la) that we have removed 
them from our sample. 



3.2 Allowed Boundaries 

As in previous studies, we require that the endpoints for 
each feature, as determined on a SN-by-SN basis by our fit- 
ting algorithm ( see Section |3.4[), be within predetermined 

D04||N 



boundaries (e.g., |Folatelli|2004| |Nordin et al.|2011b[ ). These 
boundaries allow the fitting routine to make sure that we 
are considering only a single feature at a time and not also 
including a neighbouring feature. The size of these bound- 
aries is necessary to account for variations in spectral fea- 
ture width and expansion velocity among SNe, as well as the 
temporal evolution of these values. Even though the formal 
boundaries of neighboring features (among the totality of 
SN spectra) may overlap, the actual endpoints of neighbor- 
ing features in a given spectrum will not. We also note that 
not every pixel in every spectrum will be part of a feature; 
some will effectively correspond to "continuum" flux. 

The boundaries used here are slightly modified from 
those in previous studies; our boundaries tend to be wider 



than those of our predecessors (e.g., Folatelli||2004 Nordin 
et al.||2011b| >. The final values of our boundaries were a re- 
sult of extensive testing and represent a compromise between 
including as many fits as possible (by accounting for differ- 
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ences in spectral feature width and velocity) while making 
sure that no endpoints overlap with a neighboring feature. 
The boundaries used in this work for each spectral feature 
can be found in Table [l] 



3.3 Initial Steps 

For each spectrum we begin by removing the host-galaxy re- 
cession velocity and correcting for Milky Way (MW) redden- 
ing according to the values presented in Table 1 of BSNIP I 
and assuming that the extinction f ollows the|Cardelli et ah] 
(19891 extinction law modified by lO'Donnelll (119941). The 



spectrum is then smoothed using a Savitzky-Golay smooth- 



ing filter ( Savitzky fc Golay|l964[ ). By smoothing each spec- 
trum we can effectively remove errant flux values in individ- 
ual pixels (or pairs of consecutive pixels) due to, for example, 
cosmic rays. Measured parameters from a random subset of 
spectra before and after this smoothing are effectively equal. 

From this point in the procedure we only focus on the 
spectral region near the current feature being measured. The 
S/N is then calculated and no attempt is made to measure 
the spectral feature if the S/N is less than 6.5 pixel -1 over 
the entire feature. This cutoff is based on the fact that no 
data with S/N < 6.5 pixel -1 yielded reasonable spectral fits. 
We also calculate the uncertainties in the measured flux at 
each pixel by taking the root-mean square error (RMSE) of 
wavelength bins centred on each pixel. 



3.4 The Pseudo-Continuum 

One of the most difficult aspects of a study such as this 
is defining suitable continua for the spectral features. Since 
SN la spectra consist of broad, heavily blended absorption 
features, the physical spectral continuum is nearly impossi- 
ble to define accurately ( jNordin et al.|[2011b| ). However, we 
can define a pseudo-continuum for each feature which will 
allow us to measure spectral features accurately and con- 
sistently, although the direct physical interpretation of such 
measurements is complicated and beyond the scope of this 
paper ( |Folatelli||2004[ ). 

In order to define the pseudo-continuum, the local min- 
imum of the data for the current spectral feature is deter- 
mined. The local slope of the data is then calculated in wave- 
length bins to either side of this minimum until the slope 
changes sign (i.e., we have reached a local maximum). The 
Mg II, Fe II, and S II "W" features consistently have local 
maxima within the features themselves and thus our usual 
method will determine the endpoints incorrectly. Therefore, 
we began calculating the slope of these features just inside 
the inner edges of their allowable ranges (see Table [I). Fur- 
thermore, the flux blueward of the O I triplet rarely reached 
a local maximum before entering the region surrounding the 
Si II A6355 feature. Again, this would lead to an inaccurate 
pseudo-continuum definition. To remedy this, we allowed the 
blue endpoint of the pseudo-continuum of this feature to be 
defined where the slope of the flux is > 0.002 (since it rarely 
actually changes sign, which is the endpoint criterion for all 
other spectral features). 

Once these two endpoints are determined, a quadratic 
function is fit to the data in wavelength bins centred on each 
endpoint. If either fit results in a concave upward parabola, 



we consider the endpoints to be ill-determined and no fur- 
ther attempt to fit the feature is made. In addition, if either 
parabola's peak is outside the allowed boundary range for 
the feature being measured (see Table [l]), we consider the 
pseudo-continuum to be incorrectly defined and again no 
further attempt to fit the feature is made. However, if both 
fits resulted in concave downward parabolas, with peaks 
within the allowed boundary range for the spectral feature 
in question, we connect the peaks of each parabola with a 
line and define this as the pseudo-continuum. 

This method is similar to those used by previous stud- 
Blondin et aL]|2011| |Nordin et al.||2011b[ ), though 



(e-g- 



our use of a quadratic fit to the region near each endpoint 
is somewhat unique. This extra step can be thought of as 
an additional local smoothing function to ensure that each 
pseudo-continuum endpoint truly represents a local maxi- 
mum in the flux and not simply the top of a noise spike that 
remains in the data even after our initial smoothing. Also, 
note that our pseudo-continuum definition is completely au- 
tomated (i.e., the endpoints are not manually chosen). 

When a pseudo-continuum is determined for a given 
spectral feature, we record the flux at the blue and red end- 
points of the feature (F), and F r , respectively), which are 
effectively the peaks of the two parabolas mentioned above. 
These values correspond to ftbiuc and ft. pe ak (respectively) 



from Blondin et al. (2011). The uncertainties of these values 



are simply the calculated RMSE at these pixels. An example 
of F), and F r (along with the pseudo-continuum and the rest 
of our spectral measurements) is shown in Figure [5] 

For all features with a well determined pseudo- 
continuum, we also calculate the pseudo-equivalent width 
(pEW; e.g., Garavini et al.|20'07 ) defined as 



pEW = J2 AA » 



fc(Xi) - fjXi) 
UK) 



(l) 



where Xi are the wavelengths of each pixel in the spectrum 
ranging from the blue endpoint to the red endpoint (as de- 
fined by the pseudo-continuum), AAi is the width of the i th 
pixel, /(Ai) is the data's flux at Xi, and f c (Xi) is the flux of 
the pseudo-continuum at A^. The lcr uncertainty of the pEW 
was calculated by error propagation of the uncertainty in the 
measured flux at each pixel. Somewhat surprisingly, varying 
the exact choice of pseudo-continuum endpoints added neg- 
ligibly to the uncertainty of the pEW (as well as all of the 
other values measured). The pEW is represented schemati- 
cally in the bottom panel of Figure [2] 



3.5 Spline Fits 

Once a pseudo-continuum is calculated for a given spectral 
feature, a cubic spline is fit to the smoothed data between 
the endpoints previously determined. However, no attempt 
is made to fit any of the Mg II or Fe II features in this man- 
ner. These complexes consist of so many blended spectral 
lines that it is unclear which reference wavelength to use 
when attempting to define an expansion velocity (which is 
one of the main reasons for fitting a function to the data in 
the first place). 

Other functional forms were considered before the 
spline was chosen. This included the Gauss-Hermite (Ivanl 



der Marel & Franx 19931, Gaussian, and sixth-order poly- 



BSNIP II: Initial Spectral Analysis 5 




5900 6000 6100 6200 6300 
Rest Wavelength (A) 



relativistic Doppler equation, we calculate the expansion ve- 
locity (v). Thus, all velocities used in this study are relative 
to the deepest component of each spectral feature. As men- 
tioned previously, varying the pseudo-continuum endpoints 
did not significantly change the measured values of A m i n , 
and thus we impose a 2 A uncertainty on the wavelength at 
which the spline fit reaches its minimum. We then propa- 
gate this error through the relativistic Doppler equation to 
calculate the uncertainty in the expansion velocity. 

The spectral feature being measured is then normalised 
to the pseudo-continuum, and both the relative depth of 
the feature (a) and its full width at half-maximum inten- 
sity (FWHM) are computed. The bottom panel of Figure [5] 
shows both a and the FWHM. The uncertainty of a is 
taken directly from the calculated RMSE (normalised to 
the pseudo-continuum). The uncertainty of the FWHM is 
the standard deviation of the FWHM values when varying 
the half-maximum level. 



o 

z 




5900 6000 6100 6200 6300 
Rest Wavelength (A) 

Figure 2. An example of the spectral measurements used in 
this study. The feature shown is Si II A6355 from a spectrum of 
SN 2002he at maximum brightness. In both panels the data are 
the solid black curve, the spline fit is the solid red curve, and the 
pseudo-continuum is the dashed black line. The top panel shows 
the SN spectrum along with the flux at the endpoints of the fea- 
ture (Fj, and F r ) as well as the minimum of the spline fit (A m i n ) 
which is used to calculate the expansion velocity (v). The bottom 
panel shows the spectrum normalised to the pseudo-continuum, 
in addition to the relative depth (a) and FWHM of the feature, 
and a schematic representation of the pEW. 



normal. In many cases all of these functions matched the 
data relatively well, though the spline fits matched the data 
better in more cases than the other functions. The Gauss- 
Hermite and Gaussian parameters calculated from spectra 
of astrophysical sources can be directly related to physical 
properties of the source in question (e.g 



van der Marel & 



Franx 19931. Unfortunately, such an interpretation is un- 



realistic since a true continuum is not being measured in 
SN spectra and the features that are measured are actually 
blends of many spectral lines. 

Using the wavelength at which the spline fit reaches its 
minimum (A m i n , labelled in the top panel of Figur e [2] ), along 
with the reference rest wavelength (listed in Table lib and the 



3.6 Final Inspection 

While the aforementioned automated fitting procedure is 
quite robust, we felt it was wise for each spectral feature 
in each spectrum (for which a pseudo-continuum was deter- 
mined) to be visually inspected by more than one coauthor. 
Thus, 3141 spectral features and their fits were examined by 
eye. 

About 10-40 per cent of the time (depending on the 
feature) the pseudo-continuum endpoints passed the auto- 
mated fitting criteria but did not accurately reflect the ac- 
tual edges of the spectral feature in question. This was usu- 
ally due to the measured feature being blended with a neigh- 
boring spectral feature. In these cases we simply removed 
these fits from the rest of our analysis, since if the pseudo- 
continuum was unreliable then any of the other measure- 
ments would be as well. 

Of the features determined to have accurate pseudo- 
continua, sometimes the spline fit did not accurately re- 
produce the actual minimum of the flux. This meant that 
the calculated expansion velocities, relative spectral feature 
depths, and FWHM would be inaccurate. In cases such as 
this, these spectral measurements are ignored, but Ft, F r , 
and the pEW are recorded since these values are based on 
the accuracy of the pseudo-continuum alone and are unaf- 
fected by the spline fit to the data. As mentioned previously, 
all Mg II and Fe II data fall into this category since we do 
not even attempt a spline fit for these features. About 30- 
50 per cent of the Ca II H&K, O I triplet, and Ca II near-IR 
triplet lines were found to have unreliable spline fits. Less 
than 12 per cent of each of the three Si II features and none 
of the S II "W" data had untrustworthy spline fits. 



4 RESULTS 

The dataset used in this analysis, after the aforementioned 
automated cuts were applied and the spectra were visu- 
ally inspected, contains 432 spectra of 261 SNe la. Each 
of these has a measured pseudo-continuum for at least one 
spectral feature. A summary of the number of spectra and 
objects with well-defined pseudo-continua and the number 
with "good" spline fits can be found in Table[2] As described 
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Table 2. Number of Spectra and Objects Measured per Feature 



1 cdl U.I t. 


Good 
Pseudo-Continuum 


Good 
Spline Fit 


Spectra 


SNe 


Spectra 


SNe 


Ca II H&K 


281 


191 


172 


128 


Si II A4000 


188 


137 


172 


129 


Mg II 


219 


163 








Fe II 


313 


217 








S II "W" 


240 


179 


240 


179 


Si II A5972 


204 


156 


166 


129 


Si II A6355 


366 


239 


360 


235 


O I triplet 


192 


139 


109 


84 


Ca II near-IR triplet 


301 


201 


129 


103 



region. The measured values for the Mg II feature are shown 
in Table B3. 



4.1.4 Fell 

The Fe II feature suffers from the same blending issues as 
the Mg II feature, and thus we again do not attempt to fit 
it with a spline. Another similarity with the Mg II feature is 
that after about 7 d past maximum brightness the red edge 
of the Fe II feature becomes difficult to distinguish from the 
blue edge of the S II "W." As in the case of Mg II and Si II 
A4000, we consider the continuum to be ill-defined for both 
Fe II and S II "W" when it is unclear if the two features are 
distinct. The measured values for the Fe II feature can be 
viewed in Table B4. 



in Section |3.6[ a feature in a given spectrum that has a 
"good" spline fit will have a well-defined pseudo-continuum 
by construction, though the opposite is not necessarily true. 
All measured values for each feature can be found in the 
tables in Appendix B. 



4.1 Comments on Individual Spectral Features 

4.1.1 CallH&K 

The Ca II H&K feature usually falls completely within our 
data and we are able to accurately measure it in many of our 
spectra. Sometimes the left edge of this feature is not well 
defined due to either noise at the bluest end of our data or 
the complex spectral shape. The velocity of this feature may 
be somewhat uncertain (especially at the earliest epochs) 
due to detached, high-velocity absorption that is sometimes 
observed in Ca II H&K (e.g., |Branch et al.|20 05). The mea- 
sured values for the Ca II H&K feature can be found in 
Table Bl. 



4.1.2 Si 11x4000 

The Si II A4000 feature is measurable in many of our spectra 
before and near maximum brightness. By about 7 d past 
maximum, it often weakens to the point where it becomes 
indistinguishable from the complex blend of spectral lines we 
refer to as the Mg II feature. In spectra where it is unclear 
whether Si II A4000 is a distinct feature or blended with 
Mg II, we consider the continuum to be ill-defined for both 
spectral features. The measured values for the Si II A4000 
feature are presented in Table B2. 



4.1.5 SII "W" 

After about 7 d past maximum, the S II "W" weakens signif- 
icantly and becomes blended with both Fe II (as mentioned 
above) and Si II A5972. The two broad features that make 
up the tell-tale "W" shape of this S II feature are sometimes 
so broadened (at the highest expansion velocities) that they 
are almost indistinguishable. Note that all velocities derived 
for the S II "W" are with respect to the redder of the two 
broad features (~5624 A). The measured values for the S II 
"W" feature are displayed in Table B5. 



4.1.6 Si 11X5972 

The Si II A5972 feature, as stated previously, becomes 
blended with the S II "W" after about 7 d past maximum 
brightness. It also becomes blended with the usually much 
stronger Si II A6355 feature near this epoch as well. Fur- 
thermore, this feature can sometimes be contaminated by 
Ti II absorption, especially in Ia-91bg objects (Filippenko 



et al.|1992b I. The spectral range over which we fit the Si II 
A5972 feature also includes Na I D at rest (i.e., from the 
Milky Way) and, for most of the objects presented here, 
at the redshift of the SN host galaxy. The vast majority of 
our data do not show strong Na I D absorption from either 
source, but there are a few spectra where it is detected. No 
attempt is made to correct for this absorption or interpo- 
late over it; we simply point out that our measured pEW 
of the Si II A5972 is perhaps larger than the actual value 
in a few cases due to the added absorption from Na I D. 
The measured values for the Si II A5972 feature are listed in 
Table B6. 



4.1.3 Mgll 

As mentioned previously, we did not attempt to fit any of 
the Mg II features with a spline function. This is mainly due 
to the fact that this feature is actually made up of blends of 
many iron-group element (ICE) spectral lines and thus was 
extremely complex to fit (even with something as generic 
as a spline function). In cases where a spline would have 
fit the data fairly well, we did not record its velocity since 
it is unclear which rest wavelength to use when attempting 
to define an expansion velocity for such a complex spectral 



4.1.7 Si 11X6355 

The characteristic spectral feature of SN la spectra near 
maximum brightness is the Si II A6355 trough. Unsurpris- 
ingly, this is the line for which the most "good" fits were 
obtained (see Table I2J. As we have already pointed out, 
the nearby (though usually weaker) Si II A5972 feature be- 
comes somewhat blended with Si II A6355 by about 7 d 
past maximum. However, the Si II A6355 line is usually so 
much stronger than Si II A5972 that we are still able to ob- 
tain "good" fits to Si II A6355 well after 7 d past maximum 
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brightness. The measured values for the Si II A6355 feature 
can be found in Table B7. 



4.1.8 01 Triplet 

Perhaps the most uncertain feature we attempt to fit is the 
O I triplet (notice its relatively low numbers in Table [2| . 
Even though this is not at the reddest edge of most of our 
spectra, it is in a region that is often contaminated by large- 
amplitude fringing due to the spectrograph. In addition, it 
is usually found in a part of the spectrum that is strongly 
affected by telluric absorption. Even though corrections are 
made for both the fringing and the telluric absorption (see 
BSNIP I for more information on how these correction are 
made), there often remains significant noise in this wave- 
length region. However, the O I triplet is important for us 
to investigate here, since it has been often neglected in pre- 
vious measurements of SN la spectra. For example, one of 
the largest pre-BSNIP SN la spectral dat asets had an av- 
erage wavelength coverage of 3700-7400 A (iMatheson et al 



2008). Despite the aforementioned sources of uncertainty, 



we attempt to measure this spectral feature and are suc- 
cessful for a significant number of spectra. However, we do 
caution that the measured properties of the O I triplet are 
more uncertain than the rest of the features inspected. The 
measured values for the O I triplet are given in Table B8. 



within 2 d of each other. This adds about 20 pairs to the 
test, but also increases the median relative difference be- 
tween consecutive spectra. The difference was still on the 
order of the uncertainties reported, though the difference 
in expansion velocities was often larger than our calculated 
uncertainties. This is not as enlightening as the test with 
spectral pairs separated by one day since over the course of 
two days (near maximum brightness) the expansion velocity 
of SNe la can change by almost 100 km s" 1 d _1 (see Sec- 
tion 5.1 1, which is on the order of the median of our reported 



expansion- velocity uncertainties. 

Another, similar test was conducted using the SN la 
spectra presented by Matheson et al. (20081. Due to the 



scheduling of their telescope time, their dataset consists 
mainly of well-sampled spectroscopic time series of a handful 
of SNe la, averaging more spectra per object than in BSNIP 
(Silverman et al., submitted). When applying our measure- 
ment procedure to the Matheson et al. ( 2008) data, there are 



about 40-70 pairs of spectra with "good" fits (depending on 
the feature) separated by one day. Again, the median rela- 
tive difference between pairs is larger than or about equal 
to the median uncertainty of the measurements themselves. 
This supports the idea that the measurement procedure is 
reliable and any measured spectral changes are indeed phys- 
ical. In addition, this indicates that the uncertainties calcu- 
lated by the measurement procedure are representative of 
the actual uncertainties. 



4.1.9 Ca II Near-IR Triplet 

Finally, the Ca II near-IR triplet often falls completely 
within our data which, as mentioned above, has rarely been 
the case for previous studies similar to this one. This feature 
is difficult to measure accurately due to fringing at the red- 
dest end of many of the BSNIP spectra. Like the Ca II H&K 
feature, the velocity of the Ca II near-IR triplet may be 
somewhat uncertain (especially at the earliest epochs) due 
to detached, high-velocity absorption that has been observed 
(e.g., Mazzali et al.|2005 1. The measured values for the Ca II 
near-IR triplet are compiled in Table B9. 



4.2 Self-Consistency of the Measurements 

To investigate how self-consistent and reliable the measure- 
ment procedure used in this study is, values measured from 
spectra of the same SN obtained on consecutive nights were 
compared. Since the spectra of SNe la should not evolve 
much over the course of one day, any differences in the values 
measured should come mainly from the measurement pro- 
cedure itself (in addition to the uncertainty in the spectrum 
itself). There are about 10-20 pairs of spectra with "good" 
fits (depending on the feature) separated by one day in the 
data analysed here. We find that the median relative differ- 
ence between each pair is slightly larger than or about equal 
to the median uncertainty of the measurements themselves. 
Thus, the measured change in an object's spectrum over 
the course of one day is consistent with the uncertainties we 
report for each measurement. 

This test was redone with pairs of spectra that were 
observed within 0.5 d of each other, but there were very 
few data which met this criterion and thus no useful results 
could be obtained. It was also run with pairs of spectra 



5 ANALYSIS 

5.1 Temporal Evolution of Expansion Velocities 

Much work has been done previously on studying the ex- 
pansion velocities of the ejecta of low-z SNe la as they de- 



crease with time (e.g., Barbon et al. 1990 Branch & van 



den Bergh 1993 Benetti et al. 2005 Wang et al. 20091. It 



has been claimed that there exists a population of spectro- 
scopically normal SNe la (i.e., Ia- norm) that have higher- 
than-normal expansion velocities (as determined by the Si II 
A6355 feature) near maximum brightness, and that these 
objects might have photometric peculiarities as well (Wang 



et al. 2009 Foley fc Kasen||2011 1. Wang et al. ( 2009 \ de- 



fined high-velocity (HV) objects as spectroscopically normal 
SNe Ia with spectra within 7 d of maximum brightness that 
exhibit a velocity greater than 3<r above the average at the 
epoch when the spectrum was taken. They also found that 
this 3(J cutoff was ~11,800 km s _1 at t = d. 

As will be shown below, the scatter in Si II A6355 ve- 
locity increases drastically at ages earlier than 5 d before 
maximum brightness, so this is the lower age boundary in 
the investigation of HV objects. A histogram of the Si II 
A6355 velocities for spectra within 5 d of maximum is shown 
in Figure [3] The vertical dashed line at v = 11,800 km s _1 
is the cutoff between normal and HV objects at maximum 
brigthness. 

The average velocity of spectra with velocities less than 
11,800 km s _1 and within 5 d of maxi mum is ~10,700± 
700 km s _1 , which is consistent with what Wang et al. ( 2009 1 



found The average (linear) change in velocity with time 



2 Even though most of the spectra used here were also used by 
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10 12 14 

Velocity (10 3 km s _1 ) of Si II A.6355 



5.1.1 Call 

The evolutionary trends of the two Ca II features (Figure [4] 
top row) are quite similar to each other. Both have a large 
range of velocities at t < — 5 d, reaching values as high 
as 26,000-30,000 km s _1 , which is likely due to detached, 
high- velocity absorption that has been observed in these fea- 
tures (e.g., |Branch et al.|[2005] |Mazzali et al.||2005| ). How- 
ever, these highest velocities decrease rather quickly. For 
—5 < t < 20 d the velocities of both features are quite 
constant (at ~ 12,000-16,000 km s _1 ), with only a hint of 
decreasing with time. For both of these features, the HV 
SNe have higher typical post-maximum velocities, but the 
difference is not too significant, and the ranges of veloci- 
ties spanned by the normal and HV objects are highly over- 
lapping. Similarly, Ia-91T/99aa objects have slightly higher 
than average velocities near maximum brightness and Ia- 
91bg objects tend to have normal to low velocities. 



Figure 3. A histogram of the velocities of the Si II A6355 feature 
for spectra within 5 d of maximum brightness. The vertical dashed 
line at v = 11,800 km s _1 is cutoff between normal and HV 
objects at maximum brightness. 



of all spectra with velocities less than 11,800 km s 1 and 
t > —5 d is 38 km s _1 d _1 , again consistent with the findings 



of Wang et al. (20091 



To determine if a n object should be considered HV or 



normal, the cutoff from Wang et al. ( 2009 1, 11,800 km s , is 



applied at t — 0. The average change in velocity with time is 
then used to extrapolate that cutoff value from t = —5 d to 
t = 10 d. The lower age boundary was mentioned above, 
while the upper age boundary comes from the fact that 
the velocities of HV and normal SNe begin to significantly 
overlap by about 10 d past maximum. Three SNe, all with 
velocities close to the cutoff between HV and normal, had 
one spectrum with a velocity that was above the cutoff and 
one below the cutoff. These were ultimately classified using 
the velocity of the spectrum closest to maximum brightness. 
The results of this classification can be found in the "Wang 
Type" column of Table |A1| Objects for which a "Wang 
Type" could not be determined are either spectroscopically 
peculiar (according to SNID) or have no Si II A6355 velocity 
in the range -5 < t < 10 d. Of the 140 SNe for which a 
"Wang Type" is determined, about 27 per cent are HV (for 
comparison, ~35 per cent of the 158 objects in the sample 
of |Wang et al.|2009| were found to be HV). 

The temporal evolution of the expansion velocities for 
each of the seven spectral features with measured velocities 
is shown in Figure [4] and Figure [5] The colour of each data 
point corresponds to the "Wang type" (i.e., whether it is 
a normal or high- velocity SN or undetermined), while the 
shape of each data point corresponds to its "SNID type." 
All following figures in this work employ the same colour 
and shape scheme unless otherwise noted. 



| Wang et al.| l [2b09[ ), they used data from non-BSNIP sources as 
well. In addition, the method of measuring expansion velocities 
differs in the two studies (see Section [3] for more information on 
the measurements procedure used here). 



5.1.2 St II 

As mentioned above, the velocities of the Si II A6355 feature 
(Figure |5| span a huge range of values at t < —5 d. How- 
ever, at t > —5 d, the velocities decline relatively linearly 
with time. The typical velocity near maximum brightness 
for the la-norm is ~ 11,000 km s _1 . By construction, the 
HV objects have higher velocities than the normal objects 
for —5 < t < 10 d, but this also appears to hold at earlier 
times. At later times, the velocities of HV and la-norm ob- 
jects become quite similar (thus, there are some blue squares 
below the dashed line in Figure[5]and these epochs). Most of 
the spectroscopically peculiar objects from all subclasses (Ia- 
91bg, Ia-91T, and Ia-99aa) have lower than average veloci- 
ties near maximum brightness, with the Ia-91bg SNe having 
the lowest velocities measured in this work. Note that in 
this study we excluded some of the most peculiar SNe la, 
such as SN 2002cx-like objects (e.g., |Li et al.|2003 Jha et al. 
2006) and possible super-Chandrasekhar mass objects (e.g., 
|Howell et al.|20 06 Silve rman et al.|2011| ), which show even 
lower velocities. 

The temporal evolution of the Si II A4000 feature (Fig- 
ure [4] middle left) is very similar to that of Si II A6355, 
except it has less velocity scatter at t < — 5 d. In fact, for 
~~ 10 ~ t J$ 10 d the Si II A4000 feature appears to have a 
linear decline, and the normal and HV objects remain well 
separated during those epochs. The typical velocity of the 
Si II A4000 feature for all objects matches that of Si II A6355 
for only the normal-velocity objects. Furthermore, the high- 
est velocities seen in Si II A6355 are not exhibited by Si II 
A4000, most likely due to the fact that we are unable to mea- 
sure such high velocities for Si II A4000; at these values the 
Si II A4000 feature becomes blended with the much stronger 
Ca II H&K line. 

Conversely, the Si II A5972 feature (Figure [4] middle 
right) shows a large velocity scatter at t < —5 d and a 
significant overlap between the normal and HV SNe at t > 
5 d. The typical Si II A5972 velocity, as well as the range 
spanned, match well with Si II A6355, though perhaps also 
lacking some of the highest velocities. The upturn in many 
of the velocities at t ~ 5 d is likely due to blending between 
the Si II A5972 feature and the Na I D line which can appear 
in SN la spectra near this epoch (e.g., Branch et al. |2005[ ). 
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Figure 4. The velocity versus rest-frame age relative to maximum brightness for various spectral features: (top left) 172 spectra of 128 
SNe for Ca II H&K, (top right) 129 spectra of 103 SNe for the Ca II near-IR triplet, (middle left) 172 spectra of 129 SNe for Si II 
A4000, (middle right) 166 spectra of 129 SNe for Si II A5972, (bottom left) 240 spectra of 179 SNe for the S II "W" feature, and (bottom 
right) 109 spectra of 84 SNe for the O I triplet. Blue points are high-velocity (HV) objects, red points are normal-velocity objects, and 
black points are objects for which we could not determine whether the SN was normal or high velocity (see the text for further details 
regarding how HV SNe are denned). Squares are la-norm objects, stars are Ia-91bg objects, upward-pointing triangles are Ia-91T objects, 
downward-pointing triangles are Ia-99aa objects, and circles are objects which do not have a spectroscopic subclass (see BSNIP I for 
further details regarding how these subclasses are defined). If uncertainties in the velocities are not displayed, then they are smaller than 
the size of the data points. 
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Figure 5. The velocity of Si II A6355 versus rest-frame age rel- 
ative to maximum brightness of 360 spectra of 235 SNe. Colours 
and shapes of data points are the same as in Figure [4] The 
solid line is the mean evolution of the normal-velocity SNe with 
t > —5 d, the dotted lines are the la uncertainties, and the dashed 
line is the cutoff between HV and normal-velocity objects (only 
defined for spectra with —5 < t < 10 d). If there are multi- 
ple spectra for a given object, the spectrum closest to maximum 
brightness is used to determined whether it is normal velocity 
HV. Thus, there can be blue sqaures that fall below the dashed 
line. Uncertainties in the velocities are smaller than the size of 
the data points. 



5.1.3 SII 

The temporal evolution of the S II "W" velocity (calcu- 
lated using the minimum of the redder of the two absorp- 
tion features relative to its rest wavelength) is shown in 
the bottom left of Figure [4] and is quite linear from about 
10 d before maximum brightness through 10 d after max- 
imum. At about 5 d before maximum the typical velocity 
is ~ 10,000 km s _1 , and by 5 d after maximum the typi- 
cal velocity is ~8000 km s _1 . As expected, the HV objects 
have larger velocities on average than the normal objects, 
but again there is significant overlap between the two sub- 
classes. In addition, it seems that some of the HV objects 
decrease in velocity more quickly than the normal objects 
such that they have about the same velocity as the la-norm 
at the time of maximum brightness (see Section [5.2| for more 
information regarding the change of velocities with time). 
The Ia-91T/99aa objects have relatively large velocities for 
the S II "W" feature, similar to the HV SNe. Like the ve- 
locities of the Si II A6355 feature, the Ia-91bg objects have 
many of the lowest velocities measured in this work (where 
we have excluded some of the slowest expanding, most pe- 
culiar SNe la). 



5.1.4 01 

The O I triplet (Figure [4] bottom right) behaves in a man- 
ner similar to that of the Ca II features. The velocities re- 
main constant in most of our spectra. However, the typical 
post-maximum velocity, ~11,000 km s" 1 , is lower than that 
of the Ca II features and has less scatter. The cluster of 6 



spectra at extremely low velocities all have relatively weak 
O I triplets, and since this feature is so broad, determin- 
ing the exact minimum when it is weak can be somewhat 
inaccurate. That being said, the visual checks of these fits 
reveal that they should be considered "good" fits. Like the 
Si II A6355 feature, most of the spectroscopically peculiar 
objects from all subclasses (Ia-91bg, Ia-91T, and Ia-99aa) 
have lower than average velocities, although there is much 
scatter in the velocities of the Ia-91bg SNe. 

5.1.5 Summary of Velocity Evolution 

Within a few days of maximum brightness, the typical ve- 
locities of the Ca II features are 12,000-15,000 km s" 1 , those 
of the Si II and O I features are 10,000-11,000 km s" 1 , and 
that of the S II feature is ~9000 km s -1 . These differences in 
velocities support the idea of the layered structure of SN la 
ejecta, with Ca II found predominantly in the outermost 
(i.e., fastest expanding) layers, and O I, Si II, and S II found 
mainly in the inner (i.e., slower expanding) layers. Ca II 
shows the largest velocity scatter and is likely well-mixed 
from the outermost layers into moderately deeper layers of 
the ejecta. O I, Si II, and S II all show velocity scatter simi- 
lar to each other (and smaller than that of Ca II), implying 
that they are probably not as thoroughly mixed throughout 
the ejecta. 

When comparing the velocities of most of the features 
investigated here, the typical near-maximum velocity of the 
normal SNe is in fact lower than that of the HV objects. 
However, this difference in velocities is relatively small in 
most cases, and there is significant overlap in the veloci- 
ties spanned by the normal and HV SNe. By construction, 
the normal and HV objects have significantly different Si II 
A6355 velocities, but even in the other Si II features there 
is a fair amount of overlap among the two subclasses. Fur- 
thermore, in Figure [5] there does not appear to be a strong 
distinction between the normal and HV objects, and thus 
a sharp cut to define these subclasses does not seem very 
well-motivated. 

Ia-91bg objects tend to have expansion velocities lower 
than those of spectroscopically normal SNe, with the Si II 
A6355 velocities being the most extreme case, though there 
is much overlap among the velocities for these subclasses. Ia- 
91T/99aa objects, on the other hand, are much more com- 
plicated. While there is significant scatter in the velocities 
calculated for these SNe, they are sometimes found to have 
higher than average velocities (in the Ca II and S II features) 
and sometimes smaller than average velocities (in the Si II 
A6355 and O I features). 

Even though we question the idea of two distinct veloc- 
ity populations of spectroscopically normal SNe la, it has 
been shown that normal and HV objects may have different 



intrinsic reddening or colours (Wang et al. 2009 Foley & 



Kasen 20111. The BSNIP data seem to indicate, however, 



that instead of two distinct populations with two different 
intrinsic colours, there is likely a (nearly) continuous distri- 
bution of near-maximum velocities which may be correlated 
with intrinsic colour (or reddening). This is supported by re- 
cent theoretical models and interpretations of these models 
that explain the existence of normal and HV SNe la based 
on different viewing ang les to the SNe JKasen & Plewa|2007 



|Kasen et al.|[2009l paeda et al.|[20To] 



Foley fc Kasen||2011 
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Foley et al.|2011 l. Further investigations into these models 
and the photometric differences between HV and normal- 
velocity objects will be conducted in future BSNIP studies. 

5.2 Velocity Gradients 

One way to quantify how expansion velocities of SNe la 
evolve with time is to calculate their velocity gradient. 



Benetti et al. (20051 defined the velocity gradient, 



Av /At, as the "average daily rate of decrease of the ex- 
pansion velocity" of the Si II A6355 feature. Note that the 
BSNIP sample is not the best suited for this kind of study 
since the average number of spectra per object is relatively 
low (~2; Silverman et al., submitted). Therefore, the ma- 
jority of the objects in our dataset have only a single near- 
maximum spectrum and we do not attempt to determine a 
velocity gradient for these objects. Despite this limitation, 
there are still quite a few SNe with multiple spectra in the 
sample analysed here for which we can calculate a velocity 
gradient. 

While previous studies have mainly used only post- 
maximum velocities, for our v calculations we utilise veloc- 
ities from spectra with t ^ —5 d. This is reasonable since, 
as can be seen in Figure [5j the velocities of the Si II A6355 
feature decay linearly starting at t ~ —5 d, and this allows 
us to add more SNe to our v investigation. The velocity 
gradient is calculated using a linear least-squares fit to all 
the velocities of a given SN la measured from spectra with 
—5 ^ t ^ 20 d. For objects with more than two velocities, 
the uncertainty in v is computed from the linear fit. For 
objects with only two velocities, the uncertainty in v is cal- 
culated by propagating the uncertainties in the two velocity 
measurements themselves. 

Using the velocity gradient, |Benetti et al.| (2005) found 
that their sample of 26 SNe la could be divided into three 
subclasses. The high velocity gradient (HVG) group had the 
largest velocity gradients (v > 70 km s _1 d _1 ) and consisted 
of la-norm, while the low velocity gradient (LVG) group had 
the smallest velocity gradients and included la-norm as well 
as Ia-91T/99aa. The third subclass (FAINT) had the low- 
est expansion velocities, yet moderately large velocity gra- 
dients, and consisted of subluminous SNe la (i.e., Ia-91bg) 
with Amis(B) > 1.6 mag. The same subclasses and criteria 
for membership are used in the study presented here. Note 
that t he l argest velocity gradient presented by |Benetti et al. | 
(2005} or |Hachinger et al.| ( [2006] ) is 125 km s 1 d"\ while 
the BSNIP data contain 12 SNe with v > 200 km s _1 d" 1 . 

Despite the fact that the BSNIP dataset only averages 
about two spectra per object, a velocity gradient can be cal- 
culated for 61 of the SNe la. The computed values of v (and 
their uncertainties), along with the photometric references 
which are the sources of the Amis values used to determine 
whether or not a SN is FAINT, are presented in Table [3] 
The classification (i.e., Type) of each of these objects is also 
shown in Table [3] as well as in the "Benetti Type" column 
of Table [All 

The 11 objects in Table [3] that have classifications 
marked with a "?" are uncertain since light-curve shape in- 
formation is unavailable, so their classification is based only 
on their velocity gradient. As we will show below, FAINT 
objects have similar v values to HVG objects, thus some of 
the HVG? could in fact be part of the FAINT subclass. The 



6 SNe in Table [3] that have classifications marked with a '*' 
use the MLCS2k2 A parameter ( Jha et al.|2007 l as a proxy 
for Ami 5 [f] 

Most LVG and FAINT objects are also normal velocity 
(as opposed to HV) objects and have the lowest velocities 
observed (specifically in the Si II and Ca II features). This 
confirms what was observed by Benetti et al. (20051. The 



HVG subclass contains both normal and HV objects. Many 
HVG objects have relatively high velocities at early times, 
but then evolve to have normal to somewhat low velocities 
by about 10 d past maximum brightness. This confirms what 
was seen in at least one previous study ( Pig nata et al.|2008 1 , 
but differs from what has been seen (and assumed) in other 
previous work ( jBenetti et aL]|2005| |Wang et alT||2009| |Foley| 
fc Kasen|2011 1. In these studies, HVG objects were claimed 



to have higher velocities than LVG and FAINT objects from 
before maximum through i w 10 d. While most HVG SNe 
start out at high velocities before and near maximum bright- 



ness, our data are in agreement with Pignata et al. ( 2008 1 



and show that they evolve to average (or even relatively low) 
velocities by only a few days after maximum. 

The objects in the BSNIP data with the highest veloc- 
ity gradients, which have not been seen previously, appear to 
follow this paradigm of having high velocities before maxi- 
mum through t ~ 10 d. These observations can be explained 



using the off-centre explosion models of Maeda et al. ( 2010 1 



In their models, different viewing angles will result in a wide 
range of observed velocities and velocity gradients. Further- 
more, before maximum the highest velocity objects have the 
largest velocity gradients, but the velocities of all of their 
models become quite similar by 10 d after maximum, much 
like what is observed in our data. It is also interesting to 
note that jMaeda et al. | ( |2010[ ) have at least one model that 
we would classify with v « 300-400 km s _1 d _1 . 

In Table [1] we present the averages and standard devia- 
tions of Amis and v for each subclass (with and without SNe 
having uncertain classifications due to a lack of photometric 
data), as well as the number of SNe in each subclass. HVG 
is the largest group and only a few SNe fall into the FAINT 
subclass. Comparing these numbers to those of Benetti et al. 
( 2005 ) , we find a similar number of LVG and FAINT objects, 
but significantly more HVG SNe. These differences are inter- 
esting to note, but may not have any physical significance 
since the sample used by Benetti et al. ([2005 ) contained 
only well-observed SNe la which perhaps biased the sample 
to contain more peculiar or bright or nearby objects since 
these have historically been better observed. They also used 
spectra from a variety of sources which could introduce sys- 
tematic biases into their analysis, and they measure veloci- 
ties from their first epoch until the last epoch in which the 
Si II A6355 feature is detected. The BSNIP dataset used 
here consists of an effectively random sample of SNe la with 
spectra having t < 20 d obtained using a small number of 



3 A SN with Amis = 1-1 is defined to have A = 1 | |jha et aT| 
|2007[ ), and the relationship between Amis an d A is roughly linear 
over most of the range of observed values of A (e.g., [Hicken et ah] 
|2009[ ). In addition, many of our objects with a velocity gradient 
have both A and Amis values known, and using these we are 
able to confirm which subclass an object belongs to based solely 
on v and A. 
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Table 3. Velocity Gradients and Interpolated/Extrapolated Velocities 



SN Name 


Type a 




LC Ref. 1 






v c 






VQ 


d 






«10 d 




SN 


1989M 


HVG 


Ganeshalingam 


et 


al.H2010| 291.84 (139 


49) 


13 


19 


(0 


42) 


10 


27 


;o 


98) 


SN 


1991bg^ 


FAINT 




Hicken et al.||2069[ 


131.31 (6 


44) 


10 


50 


(0 


07) 


9 


19 


(1 


1(1) 


SN 


1994D+ 


LVG 


z 


Hicken et~aTTl |20091 


33.43 (6 


77) 


10 


60 


(0 


06) 


10 


27 


'0 


09) 


SN 


1999ac 


HVG 


Ganeshalingam 


ct 


al. 


12010 


445.41 (49 


10) 


9 


90 


(0 


13) 


5 


45 


,0 


(il) 


SN 


1999cp 


LVG 


Ganeshalingam 


ct 


al. 


12010 


32.36 (15 


43) 


10 


65 


(0 


10) 


10 


32 


'0 


07) 


SN 


1999da 


FAINT 


Ganeshalingam 


ct 


al. 


12010 


131.17 (15 


56) 


11 


05 


(0 


08) 


9 


73 


;o 


14) 


SN 


1999dq 


HVG 


Ganeshalingam 


ct 


al. 


12010 


83.64 (20 


04) 


10 


92 


(0 


07) 


1(1 


08 


(1 


22) 


SN 


1999gh+ 


FAINT 




IJha et al. 


|2006l 


46.49 (10 


13) 


11 


21 


(0 


13) 


1(1 


74 





17) 


SN 


2000dk 


FAINT 


Ganeshalingam 


ct 


al. 


)2010| 108.19 (14 


04) 


11 


30 


(0 


11) 


1(1 


21 


,0 


()!)) 


SN 


2000dm 


HVG 


Ganeshalingam 


ct 


al. 


12010 


97.14 (14 


08) 


11 


28 


(0 


08) 


1(1 


31 





12) 


SN 


2000dn 


LVG 


Ganeshalingam 


ct 


al. 


12010 


49.24 (7 


95) 


10 


18 


to 


09) 


9 


69 


,0 


07) 


SN 


2001bg 


HVG* 


Ganeshalingam 


ct 


al. 


12010 


228.90 (26 


52) 


14 


56 


(0 


44) 


12 


27 


,0 


18) 


SN 


2001da 


HVG 


Ganeshalingam 


ct 


al. 


12010 


88.17 (12 


75) 


11 


50 


(0 


09) 


1(1 


62 





1(1) 


SN 


2001en 


HVG 


Ganeshalingam 


ct 


al. 


12010 


103.48 (29 


98) 


13 


31 


(0 


38) 


12 


28 


,0 


1(1) 


SN 


2001ept 


HVG 


Ganeshalingam 


ct 


al. 


[2010 


96.18 (26 


92) 


10 


45 


(0 


15) 


9 


49 


'0 


31) 


SN 


2002bo 


HVG 


Ganeshalingam 


ct 


al. 


12010 


245.14 (8 


12) 


13 


(il 


(0 


09) 


11 


16 


,0 


07) 


SN 


2002cd 


LVG 


Ganeshalingam 


ct 


al. 


12010 


17.91 (8 


33) 


14 


83 


(0 


11) 


14 


65 


'0 


07) 


SN 


2002de 


HVG 


Ganeshalingam 


ct 


al. 


pnu 


96.88 (15 


94) 


11 


89 


(0 


09) 


1(1 


92 


,0 


12) 


SN 


2002ou 


HVG? 












119.97 (14 


63) 


11 


07 


(0 


10) 


9 


87 


,0 


10) 


SN 


2002ha+ 


LVG 


Ganeshalingam 


ct 


al. 


|2010 


15.54 (15 


55) 


10 


90 


(0 


08) 


1(1 


74 


,0 


18) 


SN 


2002hd 


HVG* 


Ganeshalingam 


ct 


al. 


12010 


116.29 (22 


08) 


11 


17 


(0 


22) 


1(1 


01 


,0 


07) 


SN 


2002het 


HVG 


Ganeshalingam 


ct 


al. 


12010 


81.01 (31 


94) 


12 


57 


(0 


06) 


11 


76 


,0 


33) 


SN 


2003he 


LVG 


Ganeshalingam 


ct 


al. 


12010 


14.49 (23 


71) 


11 


39 


(0 


15) 


11 


25 


,0 


12) 


SN 


2003iv 


HVG? 












98.60 (28 


67) 


11 


42 


(0 


14) 


1(1 


43 





18) 


SN 


2004bl 


HVG? 












70.91 (9 


33) 


11 


01 


(0 


13) 


1(1 


30 


,0 


07) 


SN 


2004dt 


HVG 


Ganeshalingam 


ct 


al. 


|2010 


269.45 (8 


34) 


14 


71 


(0 


11) 


12 


01 


'0 


07) 


SN 


2004fu 


HVG? 












211.06 (27 


37) 


12 


36 


(0 


07) 


1(1 


25 


'0 


29) 


SN 


2005M 


HVG 


Ganeshalingam 


ct 


al. 


12010 


86.36 (137 


62) 


8 


77 


1 


20) 


7 


90 


,0 


19) 


SN 


2005am 


FAINT 


Ganeshalingam 


ct 


al. 


12010 


61.11 (72 


62) 


11 


74 


(0 


40) 


11 


13 


;o 


34) 


SN 


2005bc 


LVG 


Ganeshalingam 


ct 


al. 


12010 


64.63 (23 


70) 


11 


01 


(0 


13) 


1(1 


37 





15) 


SN 


2005cf 


HVG 


Ganeshalingam 


ct 


al. 


12010 


106.69 (149 


59) 


10 


14 


(0 


20) 


9 


07 


1 


74) 


SN 


2005de 


HVG 


Ganeshalingam ct al. 


12010 


70.45 (12 


73) 


10 


86 


(0 


09) 


1(1 


15 


'0 


10) 


SN 


2005cl 


LVG 


Gancshalinga m 


ct 


al. 


12010 


13.46 (20 


10) 


10 


95 


(0 


12) 


1(1 


81 


,0 


13) 


SN 


2005er 


HVG? 












228.72 (71 


21) 


10 


02 


(0 


09) 


7 


73 


;o 


67) 


SN 


2005eq 


HVG 


Ganeshalingam 


ct 


al. 


J2010 


76.47 (37 


84) 


9 


58 


(0 


08) 


8 


82 





43) 


SN 


2( 1(1.' >ki 


LVG 


Ganeshalingam 


ct 


al. 


JJUTUl 24.19 (20 


54) 


11 


27 


(0 


12) 


11 


03 


(1 


12) 


SN 


2005ms 


HVG 


Hickcn ct al. 1 2009 1 


115.00 (8 


36) 


10 


95 


(0 


09) 


9 


80 


;o 


08) 


SN 


2005na 


HVG 


Gar 


cshalingam 


ct 


al. 


12010 


288.79 (138 


17) 


1(1 


93 


(0 


10) 


8 


04 


1 


31) 


SN 


2006Nt 


HVG 




Hicken et al. 1120091 


84.84 (8 


96) 


11 


20 


(0 


06) 


1(1 


35 


'0 


11) 


SN 


2006St 


LVG 




Hicken et~aTTl |2009) 


35.46 (6 


02) 


10 


60 


(0 


07) 


1(1 


25 


;o 


09) 


SN 


2006bqt 


HVG* 


Gan 


cshalingam 


ct 


al. 


pOlol 219.70 (15 


82) 


15 


47 


(0 


20) 


13 


27 


,0 


25) 


SN 


2006bt 


HVG 


Ganeshalingam 


ct 


al. 


]20l0| 223.70 (20 


34) 


11 


04 


(0 


07) 


8 


80 


,0 


24) 


SN 


2006dm 


HVG 


Ganeshalingam ct al. 


I2010I 131.21 (23 


40) 


11 


49 


(0 


28) 


1(1 


18 





08) 


SN 


2006ej 


HVG 


Ganeshalingam ct al. 


12010 


92.24 (15 


78) 


12 


11 


(0 


07) 


11 


19 


,0 


16) 




iiiuucu 


HVG 


Ganeshalingam 


ct 


al. 


12010 


356.38 (23 


54) 


14 


50 


(0 


32) 


1(1 






1(1) 


SN 


2006ct 


LVG 


Gancshalinga m 


ct 


al. 


12010 


16.07 (23 


52) 


9 


49 


(0 


10) 


9 


32 





11) 


SN 


2006ov 


HVG? 












82.46 (23 


77) 


12 


69 


(0 


33) 


11 


86 


,0 


11) 


SN 


2006ke 


HVG? 












111.91 (22 


80) 


9 


65 


(0 


14) 


8 


53 


'0 


13) 


SN 


2006sr 


HVG 




Hicken et al.l 


2009| 


152.37 (27 


60) 


12 


03 


(0 


07) 


1(1 


51 


,0 


28) 


SN 


2007A 


LVG? 












22.10 (10 


87) 


10 


83 


(0 


12) 


1(1 


61 


;o 


07) 


SN 


2007aft 


HVG 


Ganeshalingam 


ct 


al. 


2010 


72.01 (25 


73) 


10 


91 


(0 


07) 


1(1 


19 


,0 


27) 


SN 


2007cot 


LVG 


Ganeshalingam 


ct 


al. 


2010 


41.96 (10 


01) 


11 


43 


(0 


00) 


11 


01 


'0 


12) 


SN 


2007fb 


LVG? 












52.78 (10 


90) 


11 


46 


(0 


11) 


1(1 


93 


,0 


07) 


SN 


2007gi 


HVG? 












197.68 (20 


12) 


15 


66 


(0 


09) 


13 


69 


;o 


15) 


SN 


2007gk 


HVG? 












138.74 (6 


50) 


13 


83 


(0 


09) 


12 


44 


,0 


07) 


SN 


2007hj 


FAINT 


Ganeshalingam 


ct 


al. 


2010 


133.27 (10 


06) 


12 


26 


(0 


09) 


1(1 


92 


,0 


08) 


SN 


20071et 


HVG 


Ganeshalingam 


ct 


al. 


2010 


93.35 (12 


65) 


12 


78 


(0 


18) 


11 


84 





22) 


SN 


2008dx 


FAINT* 


Ganeshalingam 


ct 


al. 


12010 


97.91 (28 


21) 


9 


62 


(0 


IS) 


8 


64 


,0 


10) 


SN 


2008ec+ 


LVG 


Ganeshalingam 


ct 


al. 


12010 


68.61 (10 


82) 


11 


01 


(0 


09) 


1(1 


33 





14) 


SN 


2008ei 


HVG* 


Ganeshalingam 


ct 


al. 


|201o[ 114.40 (23 


97) 


15 


72 


(0 


16) 


14 


57 


,0 


11) 


SN 


2008s5° 


LVG* 


Ganeshalingam 


ct 


al. 


J2010I 11.86 (17 


85) 


9 


09 


(0 


11) 


8 


97 


;o 


11) 



Uncertainties arc given in parentheses. 

t This object has more than two near-maximum spectra that were used to calculate v. 

a Classification based on the velocity gradient of the Si II A6355 line jBenetti et al.|2005| . "HVG" — high 

velocity gradient; "LVG" — low velocity gradient; "FAINT" — faint/undcrluminous. Classifications marked 

with a "?" are uncertain since light-curve shape information is unavailable. Classifications marked with a 

"*" use the MLCS2k2 A parameter |Jha et al.|2007j as a proxy for Am^. 

b Sourcc of Am 15 (or MLCS2k2 A) value. 

c The velocity gradient is in units of km s" 1 day" 1 . 

^The velocity is calculated from the Si II A6355 line and is in units of 1000 km s~~ 1 . 
e Also known as SNF20080909-030. 
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Table 4. Summary of Velocity Gradient Subtypes 



Type Amis (mag) v (km s 1 d 1 ) vq (10 3 km s 1 ) viq (10 3 km s 1 ) # of Objects 



LVG 1.16 (0.24) 31.37 (18.91) 10.96 (1.30) 10.64 (1.32) 14 

LVG + LVG? ■•• 32.13 (18.59) 10.98 (1.22) 10.66 (1.23) 16 

HVG 1.23 (0.22) 156.12 (98.23) 11.98 (1.78) 10.42 (1.75) 29 

HVG + HVG? ■•• 152.30 (89.92) 11.98 (1.78) 10.46 (1.75) 38 

FAINT 1.80 (0.13) 101.35 (35.34) 11.10 (0.85) 10.08 (0.93) 7 



Average values are shown and standard deviations are given in parentheses. 

Average Amis values are undefined for the rows that include objects with no light-curve shape information. 
Note that the average Amis values do not include the 6 objects that use A as a proxy for Amis. 



instruments and reduced by only a few people (see BSNIP I 
for more on the homogeneity of the BSNIP dataset). 

The average Amis is about the same for all of the sub- 
classes except FAINT. This is partially by construction since 
all SNe with Amis > 1-6 mag are considered FAINT. How- 
ever, the fact that LVG and HVG objects all have effectively 
the same average Amis value implies that v and Amis are 
not correlated. This has been seen before, and previous sim- 
ilar studies found nearly identical average Amis values for 
each subclass (e.g., Benetti et al.(2005 l. 

By construction, the velocity gradients of LVG objects 
are significantly lower than those of HVG objects. Perhaps 
somewhat surprisingly, HVG and FAINT SNe have similar 
values of v. The average value of v for each subclass is effec- 
tively unchanged when objects with uncertain classifications 
are included, implying that the LVG? objects are almost cer- 
tainly all bona fide LVG SNe. However, it is unclear whether 
the HVG? are truly HVG or if they are actually part of the 
FAINT subclass. While the average v for the LVG and HVG 



objects is nearly equal to those from Benetti et al. (20051, 



the average velocity gradient for the FAINT SNe is some- 
what larger in our sample (although they are within one 
standard deviation of each other). Also, the average veloc- 
ity gradient of the HVG objects is larger than that of | Benetti] 
et al. (20051, due to the handful of SNe with significantly 



larger v values than what has been seen previously. 

As stated earlier, the BSNIP sample is not the best 
suited for this kind of study since the average number of 
spectra per object is relatively low. If we restrict ourselves 
to only objects with more than two near-maximum spectra 
(i.e., SNe marked with a in Table[3]), we are left with 6 
HVG objects, 5 LVG objects, and 2 FAINT objects (which 
used A as a proxy for Amis). One might worry, based on 
these numbers, that many of our HVG objects are a result of 
the large uncertainty introduced when calculating the veloc- 
ity gradient using only two data points. However, the aver- 
age Amis and v for each subclass is consistent with that of 
the entire sample when only using objects with more than 
two spectra. It should be noted that the average velocity 
gradient for the HVG subclass does decrease when applying 
this cut and actually becomes smaller than the one calcu- 
lated by|Benetti et al.| (|2005[). Similarly, the FAINT sub- 



class' average v becomes approximately equal to the one in 



Benetti et al. (2005) when only using SNe with more than 



two spectra. One caveat is that the fact that there are more 
than two spectra of these objects in the BSNIP dataset may 
imply that they are particularly interesting objects that are 



peculiar, intrinsically bright, nearby, or well-separated from 
their host galaxy. While these may be true and could lead 
to a bias in this subsample, all of these objects are la-norm 
except for one Ia-91bg and one Ia-99aa. 



5.3 Interpolated/Extrapolated Velocities 

Once a velocity gradient is calculated for a SN la, one can in- 
terpolate/extrapolate that gradient to determine the expan- 
sion velocity at a specified epoch. |Benetti et al. ( 2005 1 de- 
fined vio as the expansion velocity of Si II A6355 at 10 d past 



maximum brightness. Similarly, Hachinger et al. ( 2006 1 in- 
terpolate/extrapolate their expansion velocities to the time 
of maximum brightness (i.e., t — d), and so vo is defined 
here as the expansion velocity of Si II A6355 at maximum 
brightness. For each SN where v is calculated, vo and vio are 
also calculated. The uncertainties of these two velocities are 
computed by propagating the uncertainties in the linear fit 
(when more than two spectra are used) or, when only two 
spectra are used to determine the velocity gradient, by prop- 
agating the uncertainties in the two velocity measurements 
themselves. The computed values of vo and Via (and their 
uncertainties) are presented in Table[3] and the averages and 
standard deviations of these velocities for each subclass are 
displayed in Table [4] 

This interpolation/extrapolation calculation allows a 
more self-consistent comparison of the expansion veloci- 
ties of different objects. By the simple fact that we calcu- 
late nonzero velocity gradients, the expansion velocities are 
changing with time and not all of our objects were observed 
at exactly the same epochs. This procedure also enables us 
to make more quantitative statements regarding the differ- 
ences in expansion velocities among various subclasses of 
SNe la. 

As expected, at maximum brightness all objects deter- 
mined to be HV do in fact have vo > 11, 800 km s _1 and the 
opposite is true for objects determined to have normal ve- 
locities. This can be seen in the top panel of Figure|6j where 
we plot do versus v for all SNe having a measured velocity 
gradient. All blue points (HV SNe) are to the right of the ver- 
tical line at vo = 11, 800 km s _1 and all red points (normal 
velocity SNe) are to the left of it. This sanity check is encour- 
aging and implies that our HV determination (as outlined 
in Section 5.1 1 is relatively robust at maximum brightness. 



While there does not appear to be a strong correlation 
between v and vo, it does at least seem that, on average, 
objects with larger velocity gradients tend to have larger ve- 
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Figure 6. v$ versus v (top) and i>io versus v (bottom) for all 
SNe where we calculate a velocity gradient. Blue points are HV 
objects, red points are normal-velocity objects, and black points 
are objects for which we could not determine whether the SN was 
normal or high velocity. Squares are low velocity gradient (LVG) 
objects, circles are possible LVG objects, stars are FAINT ob- 
jects, upward-pointing triangles are high velocity gradient (HVG) 
objects, and downward-pointing triangles are possible HVG ob- 
jects. The vertical dashed line in the top panel is our HV cutoff 
(vq = 11,800 km s -1 ). The solid line is a linear fit to the "mod- 
erate decliners" (1 < Amis < 1-5 rrmg, whose data points are 
circled). The vertical dashed line in the bottom panel is our HV 
cutoff at t = d (11,8000 km s^ 1 ) decreased by the average ve- 
locity gradient for the normal- velocity SNe (38 km s _1 d" 1 ) for 
10 d (i.e., via = 11, 420 km s _1 = 11, 800 - 10 X 38). The slanted 
dotted line in the bottom panel is our HV cutoff at t = 10 d as a 
function of v (i.e., via = 11, 8000 — 10 X ii). 



locities at maximum light, though there are plenty of SNe la 
that do not follow this correlation. The average vo for LVG 
and FAINT SNe are approximately equal, while they are 
both lower than the average vo for HVG SNe (see the fourth 
column of Table [4} . 

The solid line in Figure [6] is a linear least-squares fit to 
the "moderate decliners" (i.e., 1 < Amis < 1-5 mag, whose 



data points are circled in the Figure) of the form 
v — avo + p. 



(2) 



We calculate a = 55.0 ± 1.8 and /3 = -539 ± 22 for v in 
units of 10 3 km s _1 and v in units of km s _1 



Foley 



et al. (2011 1 fit a similar relationship to their data in order 



to derive a family of functions. This allows them to calculate 
a velocity of the Si II A6355 feature at maximum brightness 
given a spectral age and a velocity at that age. They find 
a — 32.2 and (3 — —285, which are both significantly differ- 
ent than what is calculated for the BSNIP data. The large 
amount of scatter around the solid line in Figure [6] casts 
doubt on how useful the family of functions proposed by 



Foley et al. (20111 can actually be in calculating velocities 



at maximum brightness. 

Somewhat surprising, however, is the wide range of vo 
values spanned by each of the subclasses. Even at maximum 
brightness, where the difference between HV and normal 
objects is defined, a HVG SN may have normal (or even rel- 
atively low) velocity. This is contrary to many previous stud- 
ies that often assume a one-to-one correlation between HV 
and HVG and, similarly, between normal velocities and LVG 



(e.g., Haching er et al.(200 6 Pigna taet al.|20 08 Wang et al 



2009). Kolmogorov-Smirnov tests were performed on the vq 



values of LVG and HVG SNe (both including and excluding 
objects with uncertain classifications), and we find that they 
likely come from different parent populations (p ~ 0.03). 
Thus, it is still reasonable to associate LVG SNe with nor- 
mal velocities at maximum and HVG SNe with HV objects, 
but we caution that this may not be as robust an association 
as was previously thought. 

The connection between HVG and HV is even more ten- 
uous by 10 d after maximum brightness. In the bottom panel 
of Figure|6]we plot vio versus v for all SNe where we calculate 
a velocity gradient. The vertical line is our HV cutoff value 
at t = d (11,8000 km s _1 ) decreased by the average ve- 
locity gradient for the normal-velocity SNe (38 km s _1 d _1 ) 
for 10 d (i.e., v w = 11,420 km s _1 = 11,800 - 10 x 38). 
Naively, if one measured vio < 11,420 km s _1 for a given 
SN, they might classify it as a normal-velocity object. How- 
ever, about half of our HV SNe fall in this regime. This 
is actually expected since our HV definition only included 
velocities within 5 d of maximum, so extrapolating this anal- 
ysis to 10 d past maximum may not be valid. If we instead 
plot our HV cutoff at t = 10 d as a function of v (i.e., 
v\o = 11, 8000 — 10 x v), then we get the slanted line in the 
bottom panel of Figure [6] Now, the HV SNe all fall above 
this line while the normal-velocity objects are all below it. 
Effectively, some of the HVG objects have decreased their 
velocity fast enough to "catch up" with the velocities of the 
normal-velocity objects by this epoch. By 10 d past maxi- 
mum brightness, a single velocity measurement alone is not 
sufficient to determine whether an object should be consid- 
ered HV or not. 

All three velocity gradient subclasses span nearly the 
full range of vio values and the average vio is effectively the 
same for all subclasses (see the fifth column of Table B. 
Kolmogorov-Smirnov tests were performed on the Wio val- 
ues of LVG and HVG (both including and excluding objects 
with uncertain classifications), and we find no evidence that 
they come from different parent populations. Thus, by 10 d 
past maximum brightness the distributions of expansion ve- 
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locities among LVG and HVG objects are consistent with 
each other. 

As mentioned in Section |5.2[ the off-centre explosion 



models of Maeda et al. ( 2010 1 may naturally explain the ex- 



istence of HVG SNe. They also show that models with the 
largest velocity gradients have the highest velocities near 
maximum brightness, but by about 10 d after maximum 
the expansion velocities of almost all of their models be- 
come quite similar (independent of initial velocities or ve- 
locity gradients). These models seem to have observational 
grounding in the data we present here. Further comparisons 
to the models and predictions of |Maeda et al.| ( |2010[ > , es- 
pecially at later epochs, will be made in future BSNIP pa- 
pers. Also, the velocity gradients, classifications, and inter- 
polated/extrapolated velocities discussed here will be com- 
pared to photometric properties (such as light-curve shape 
and decline rate) in BSNIP III. 



5.4 Temporal Evolution of Pseudo-Equivalent 
Widths 

As with expansion velocities of SNe la, much work has been 
done previously on studying the pseudo-equivalent widths 
of various spectral features as they change with time (e.g., 



Folatclli 2004, Garavini et al.| 


2007| |Bronder et al. 


2008 


Walker et al.|2011 


Nordin et al. 


|2011b||Konishi et al.| 


20111 



The temporal evolution of the pEW for each of the nine 
spectral features is shown in Figure [7] and Figure [8] The 
"SNID type" is displayed rather than the "Benetti type" 



(see Section 5.2 1 in the pEW figures because the velocity 
gradient of an object does not appear to be well correlated 
with its pEW measurements. 

Also shown in Figure [7] and Figure [8] are fits to the 
pEW evolution for each spectral feature (solid line) along 
with the standard error of the fit (grey region), using only 
la-norm (with normal velocities). For features whose pEW 
appears to evolve linearly with time a linear function is fit 
to the data, while features whose pEW have a sharp change 
in their temporal evolution have a quadratic function fit to 
them. This differs from previous studies which have modeled 
the behaviour of pEWs with time either using cubic splines 



(e.g., Garavini et al. 20071 or logistic functions instead of 
quadratic ones ( |Nordin et al.|2011b| ). As seen in the figures, 
the pEW temporal evolution for each spectral feature can 
be fit relatively well with either a linear or quadratic func- 
tion. Experiments with fitting the pEW evolution with cubic 
splines and logistic functions were carried out, but the fits 
were either worse or comparable to those of the linear and 
quadratic functions (which have fewer free parameters). 



Following Nordin et al. (2011b I, we use our fits of the 



temporal evolution of the pEW for each spectral feature to 
attempt to "remove" the age dependence of the pEW. To 
do this, an epoch-independent quantity called the "pEW 
difference" (ApEW) is defined; it is simply the measured 
pEW minus the expected pEW at the same epoch using 
the linear or quadratic fit. The uncertainty in ApEW comes 
from combining the uncertainty of the pEW measurement 
with the standard error of the fit. The ApEW values and 
their uncertainties can be found in Tables B1-B9. Note that 
while the fits were defined using only la-norm, ApEW values 
are calculated for SNe of all spectral types. 



5.4.1 Call 

The Ca II H&K feature and the Ca II near-IR triplet show a 
cluster of spectra with relatively large pEW at t < —5 d in 
Figure [7] (top row). This is perhaps due to detached, high- 
velocity absorption blending with the normal-velocity com- 
ponent (e.g., |Branch et al.|[2005l>. For t > -5 d, the pEW 



of the Ca II H&K feature decreases slightly and the scatter 
in the pEW values decreases markedly near t > 10 d. This 
feature is fit with a quadratic in order to encompass the 
relatively large number of objects with high pEW values at 
the earliest times as well as the decrease and eventual flat- 
tening out at later times. It is interesting to note that the 
typical pEW for the HV objects is larger than that of the 
normal- velocity objects (for t < 10 d). 

Despite the fact that there are not that many spec- 
troscopically peculiar objects plotted in the top-left panel 
of Figure [7] it seems that Ia-91bg follow the evolution of 
the la-norm objects while Ia-91T/99aa are below the typi- 
cal pEW values (and thus have negative values of ApEW). 
This matches what has been observed previously in other 



low-redshift datasets (Garavini et al. 2007 Bronder et al. 
|2008| ). 

On the other hand, the Ca II near-IR triplet pEW values 
increase linearly with time after about 5 d before maximum 
(which is why the evolution of this feature is fit with a linear 
function). Also distinct from the Ca II H&K feature, the 
typical pEW of the HV and normal objects is about the 
same. While the Ia-91T/99aa objects are certainly below 
the normal pEW values, the Ia-91bg objects have pEWs 
that are well above the normal evolution. Once again, this 
has been noted previously, though the few large pEW values 
at early times have not been seen and our sample has far 
more data points than earlier work (Folatclli 2004). 



5.4.2 St II 

The pEW of the Si II A6355 feature (Figure [FJ middle left) 
linearly increases with time before t « 10 d and shows a hint 
of a sharp upturn thereafter (however, this is likely due to 
the feature becoming blended with Si II A5972 at these later 
epochs). Also, before ~10 d past maximum the HV objects 
have a larger typical pEW than the normal- velocity objects. 
Like the Ca II H&K feature, the Ia-91T/99aa objects fall 
well below the normal evolution while the Ia-91bg span a 
large range of pEW values (from well below to well above the 
average evolution). This behaviour is similar to that seen in 
the low-z (and moderate- z) samples of Folatelli ( 2004 1 and 
Konishi et aL] poTTj ) 



The temporal evolution of the pEW of the Si II A4000 
feature (Figure [7J middle right) is quite unique. There is 
evidence for two distinct evolutionary tracks: one rising until 
2-3 d past maximum and then declining, and one constant 
(and lower) until 2-3 d past maximum and then rising. The 
two groups are effectively blended into one another by ~5 d 
past maximum. This creates a gap at relatively low values 
of pEW from a few days before maximum until a few days 
after maximum. The two-component evolution has been seen 
in other datasets, but the gap in the present sample is not 
nearly as pronounced as in some of the earlier studies which 
used many fewer data points ( Folatclli 2004J | Bronder et al.| 
|2008[ ). 
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Figure 7. The pEW versus rest-frame age relative to maximum brightness for various spectral features: (top left) 281 spectra of 191 
SNe for Ca II H&K, (top right) 301 spectra of 201 SNe for the Ca II near-IR triplet, (middle left) 366 spectra of 239 SNe for Si II A6355, 
(middle right) 188 spectra of 137 SNe for Si II A4000, (bottom left) 204 spectra of 156 SNe for the Si II A5972 feature, and (bottom right) 
219 spectra of 163 SNe for the Mg II complex. Colours and shapes of data points are the same as in Figure[4] The solid curve is a linear 
or quadratic fit to the data using only la-norm; the grey region is the standard error of the fit. 
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Figure 8. The pEW versus rest-frame age relative to maximum brightness for various spectral features: (top left) 313 spectra of 217 
SNe for the Fe II complex, (top right) 240 spectra of 179 SNe for the S II "W" feature, and (bottom left) 192 spectra of 139 SNe for the 
O I triplet. Colours and shapes of data points are the same as in Figure [4] The solid curve is a linear or quadratic fit to the data using 
only la-norm; the grey region is the standard error of the fit. 



There is no significant difference in pEW of the Si II 
A4000 feature between HV and normal-velocity objects, 
though almost all of the HV objects are found in the "upper 
evolutionary track." While there are only a small number of 
Ia-91bg objects for which we measure a pEW for the Si II 
A4000 feature, they also fall well within the "upper evolu- 
tionary track." The Ia-91T/99aa SNe are found only in the 
"lower evolutionary track." Note that la-norm objects are 
found in both tracks. 

The ApEW values of the Si II A4000 have been found to 



correlate with SN colour as well as velocity gradient ( Nordin 
|et al.||2011a| ). In BSNIP III the relationship between both 
ApEW and pEW and SN colour for this feature will be in- 
vestigated. In Figure [9] we plot ApEW of Si II A4000 against 
—i) (the minus sign is used here in order to match the veloc- 
ity gradient definition of |Nordin et al.|2011a \ for all objects 
with ^ t ^ 8 d. Our plot contains 31 SNe la as com- 
pared to the 20 objects shown in Figure 3 of |Nordin et al.| 
(2011a I, and we also follow their convention of taking the 



mean ApEW value for objects with multiple spectra in the 
epoch range studied. 

The basic trends seen in Figure [9] are unchanged if we 
use pEW instead of ApEW, and are similar to what was 
observed by Nordin et al. < 2011a[). The bi ggest difference 
between the two studies is that Nordin et al.| ( [20~lla[ ) use only 
spectroscopically normal SNe (thus they would not have the 
four black points at the bottom right of Figure |9| , and they 
do not see the most extreme HVG objects which appear 
in our dataset (i.e., the two left-most points in Figure [9|. 
Nordin et al. (2011a I claim a "strong correlation" between 



ApEW values of Si II A4000 and velocity gradient (quoting 
a Spearman rank coefficient of —0.73). A fit to the BSNIP 
data (excluding the 6 aforementioned outlier objects which 
are ignored or not seen in their sample) yields a Spearman 
rank coefficient of —0.54, which implies that the supposed 
correlation may not actually be all that significant. 

There is a relatively large scatter in the pEW values 
measured for the Si II A5972 feature (Figure [7] bottom left). 
Evidence suggests that the pEW values are trending slightly 
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Figure 9. The ApEW of Si II A4000 versus —v (the minus sign 
is used here in order to match the velocity-gradient definition of 
|Nordin et al.|2011a[ l for all 31 objects with < t < 8 d. Colours 
and shapes of data points are the same as in Figure [1] 



upward with time; however, this is driven mainly by points 
at t > 5 d where we might expect the Si II A5972 feature 
to start blending with the Na I D line, which can appear in 
SN la spectra near this epoch (e.g., Branch et al|2005 l. Ig- 
noring points at t > 5 d, the temporal evolution is relatively 
constant, although the HV objects are perhaps decreasing 
slightly with time. Once again, the Ia-91T/99aa objects have 
relatively small pEW values while the Ia-91bg SNe lie well 
above the average evolution. This matches the trends seen 



in the low-z data presented by Folatelli ( 2004 1 and Konishi 
[etTaLlpoTl] . 

5.4.3 Mgll 

The temporal evolution of the pEW of the Mg II complex 
(Figure [7J bottom right) has relatively small scatter and 
linearly increases slightly with time. The HV objects have 
larger pEW values (and more scatter) than the normal- 
velocity objects. Interestingly, this evolution is markedly dif- 
ferent from what has been seen in some other low-z samples 
( Folatelli|2004||Garavini et al|2007) . While the BSNIP data 
match these previous studies until t w 10 d, there is no ev- 
idence for the sudden increase in pEW values for this fea- 
ture. In fact, no attempt is made to measure the pEW of 
the Mg II complex beyond t « 10 d because it becomes 
too blended with the Si II A4000 feature (as pointed out by 
Garavini et al.|2007 |. Folatelli (2004) simply define a larger 
wavelength range at these epochs (see their Fig. 1, feature 
3), which will certainly increase the measured pEW. Our 
temporal evolution does, however, match what was seen by 



Walker et al. (20111 



The Ia-91T/99aa objects are yet again found to have 
lower than average pEW values, though not too low. On the 
other hand, the Ia-91bg objects have somewhat higher than 
average pEW values, but there are very few of these types 
of objects for which a pEW is successfully measured. This is 
due to the fact that when the Mg II complex is strong before 
and near maximum brightness, as it often is for Ia-91bg SNe 



thus no attempt is made to measure its pEW. This is further 
supported by the fact that we s ee no objects w i th a pEW 
of Mg II > 140 A, while data in iGaravini et all d2007h and 



the low-z sample presented by Bronder et al. ( 2008 1 contain 



several objects with such large pEW values. 



5.44 Fe II 



While the Mg II complex has a fairly tight linear pEW evo- 
lution, the pEW of the Fe II complex (Figure |8j top left) 
shows an extremely tight temporal evolution which is al- 
most certainly nonlinear (thus we fit it with a quadratic 
function). Another difference between these two features is 
that the Mg II complex pEW values are nearly constant for 
— 10 i: t IS 10 d (increasing by < 20 A over that range) 
while the Fe II complex pEWs increase dramatically during 
the same span of time (~130 A). Despite these differences, 
the data for both the Mg II and Fe II complexes fall al- 
most completely within the standard error of their respec- 
tive fits. Furthermore, the HV objects in both features have 
larger pEW values, as well as more scatter, than the normal- 
velocity objects. 

As with the pEWs of the Mg II feature, the Ia-91T/99aa 
objects are all below the fit to the la-norm. It is even more 
apparent for the Fe II complex than it was for the Mg II 
complex that the Ia-91bg objects are all found above the fit 
to the data. These differences, as well as the overall trends, 
have been recognised in numerous previous studies of both 
low-z and moderate-z SNe la (e.g., Folatelli|2004 Garavini 



et al|2007||Nordin et al~1|2011b| [Konishi et al.||2011[ ). 



5.4.5 SII 

The temporal evolution of the pEW values for the S II "W" 
feature (Figure [8] top right) is unique. The pEW of the 
majority of objects increases until maximum brightness and 
then decreases in a nearly symmetric way. A quadratic func- 
tion centred near t — d fits the data (especially the la-norm 
objects) very well. There is a possibility that the HV objects 
have higher pEW values before maximum (with equal pEW 
values after maximum), but the significance of this differ- 
ence is relatively low. The Ia-91T/99aa objects have even 
lower pEW values than the normal-velocity objects before 
maximum, though they also seem to evolve to more average 
values after maximum brightness. The Ia-91bg objects, on 
the other hand, are almost all significantly below the av- 
erage evolution at all epochs (though there are a few with 
relatively normal pEWs). The temporal evolution seen here 



due to additional absorption from Ti II (e.g., Bronder et al 



confirms what was found by Nordin et al. (201 lb I and Kon- 
|ishi et ah] ( |2011| ) in their low-z and moderate-z data. 



5.4.6 OI 

The O I triplet's average pEW values evolve in much the 
same way as those of the S II "W," with the normal pEW 
values rising until ~4 d past maximum and then declining 
(Figure [8] bottom left). There is a hint that the HV objects 
might have smaller pEW values than the normal-velocity 
objects, though the scatter in both of these subclasses is 
quite large. As mentioned in Section |4.1.8| this feature is 
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often strongly affected by telluric absorption, and thus some 
of the points which lie significantly away from the average 
evolution could be explained by imperfect telluric absorption 
corrections. This could also be the reason why we are able 
to measure pEWs of the O I triplet for relatively few Ia- 
91T/99aa objects. As was seen in the Ca II features, the 
Ia-91T/99aa SNe have relatively low pEW values, and if 
their O I triplet pEWs are also small, then the feature could 
get "lost in the noise" after an inaccurate telluric absorption 
correction. The Ia-91bg objects, however, have large pEW 
values for ( < 2 d, but quickly decrease to more average 
values at later epochs. 



5.4-7 Summary of pEW Evolution 

As mentioned above, we chose to present the "SNID type" 
rather than the "Benetti type" in the pEW versus time fig- 
ures because the velocity gradient of an object is not well 
correlated with any pEW measurements. The one species 
with a possible correlation is Si II. In all three of its features, 
while the LVG and HVG objects show very little difference 
in pEW values, there is an indication that the pEWs of the 
FAINT objects are larger than the average values at each 
epoch. This is unsurprising since the Ia-91bg objects also 
have larger-than-average pEW values for the Si II A4000 
and Si II A5972 features, and these SNe are also often un- 
derluminous (which is exactly how the FAINT subclass is 
defined). 

When comparing the pEW values of normal and HV 
objects, for most of the features investigated here the typical 
pEWs of the HV SNe are larger than those of the normal- 
velocity objects. The HV SNe also tend to have more scatter 
in their pEW values. These two points are most noticeable 
in the Ca II H&K, Mg II, Fe II, and Si II A6355 features. 
However, some of the features do show very similar pEW 
values between HV and normal- velocity objects, and the O I 
triplet's typical pEW for HV SNe is perhaps smaller than 
that of normal-velocity objects. By t ~ 10 d the average 
pEW values of HV and normal-velocity SNe are nearly equal 
in almost all of the spectral features measured. 

The velocity of the Si II A6355 feature is used to dif- 
ferentiate between HV and normal objects, so it may be 
unsurprising that the pEW of this same feature shows a 
marked split between these two subclasses. Somewhat more 
interesting is the fact that the other three features which ex- 
hibit a strong split in pEW between HV and normal- velocity 
objects are three of the four bluest features we measure 
(Ca II H&K, Mg II, and Fe II)QHV SNe may have differ- 
ent intrinsic colours which could be caused by line blanket- 
ing due to their relatively large velocities, especially at the 



alter the effective optical depth at these wavelengths, espe- 
cially for IGEs such as Fe and Mg ( Dommguez et aL||2001 



shortest optical wavelengths ( |Foley fe Kasen 2011[ ) . This is 
exactly where the greatest pEW difference between HV and 
normal-velocity SNe is observed. An alternative, yet related, 
explanation for increased pEWs for HV objects at the blue 
end of the optical regime is that a change in metallicity can 



4 The second bluest feature we measure, Si II A4000, is relatively 
weak and has a large scatter in pEW values (perhaps even show- 
ing evidence for two separate evolutionary tracks). Thus, it is 
reasonable that the difference between HV and normal-velocity 
objects is difficult to observe in the pEW values of this feature. 



Timmes et al.|2003| |. This may hint at a possible correlation 
between expansion velocity and pEW (at least for the four 
spectral features mentioned), and we investigate this possi- 
bility further in Section [5.8| Correlations between the pEWs 
of HV and normal-velocity SNe and their intrinsic colours 
(and reddening) as well as host-galaxy metallicities will be 
investigated in future BSNIP studies. 

Ia-91T/99aa objects consistently have the lowest pEW 
(and ApEW) values at all epochs of any of the SN la spec- 
tral subtypes, with the exception of the S II "W" where 
they are the lowest for t < d but increase to more average 
values after maximum. This generic trend of Ia-91T/99aa 
objects having small pEWs has been seen before (e.g 



Fo- 



|latelli||2004] |Bronder et~al~1|2008| ). The relative weakness of 
the spectral features in these SNe can be attributed to the 
fact that they tend to be overluminous and have higher in- 
ferred temperatures (e.g., Nugent et al.|[l995 l. 

On the other hand, the Ia-91bg SNe usually have the 
largest pEWs (and ApEWs) at all epochs (though the S II 
"W" feature is again the exception in that the Ia-91bg SNe 
have below average pEW values at all epochs). This basic 
trend is most readily explained (partially) by a temperature 
effect since Ia-91bg objects are usually underluminous and 
often found to be cooler than la-norm and Ia-91T/99aa SNe 
(e.g., Nugent et al. 1995 1. As mentioned above, there is often 



evidence for Ti II absorption in Ia-91bg objects which will 
blend with what has been defined here as the Mg II complex. 
However, when Mg II is observed to be this strong, it is often 
also blended with the Si II A4000 feature and a pEW would 
not be measured for that spectrum. This explains the low 
number of Ia-91bg objects for which a pEW is measured 
for Mg II. Comparisons of spectroscopic subtypes and their 
pEW values to light-curve shapes and luminosities will be 
undertaken in BSNIP III. 



5.5 Spectral Classification Using 
Pseudo-Equivalent Widths 

The differences among the various spectroscopic subtypes 
can be more quantitatively investigated by directly com- 
paring pEW values of different spectral features. |Branch| 
et al. (2006) presented a means of spectroscopically classi- 
fying SNe la using near-maximum (— 3 ^ f ^ 3 d) pEWs of 
Si II A6355 and Si II A5972 (W(6100) and W(5750) in their 
notation, respectively). This original sample has since been 
updated by Branch et al. (20091. Based on their EW mea- 



surements they split their sample into four distinct groups: 
core normal (CN), broad line (BL), cool (CL), and shallow 
silicon (SS). However, they point out that the SNe seem to 
have a continuous distribution of pEW values, and so how 
the exact boundaries are defined is not critical. 

The pEW values we measure for Si II A6355 and Si II 
A5972 for —5 t ^ 5 d, along with the boundaries from 
Branch et al. ( 2009[ ) and the median pEW value for each 
feature, are shown in Figure [To] The median uncertainty in 
both directions is shown in the upper-right corner of the 
figure. Our plot of 89 SNe la covers a parameter space sim- 
ilar to that of the 59 SNe la plotted in Figure 2 of | Branch] 
|et al.| ( |2009[ ), and ten of the eleven SNe in both datasets 
are classified as the same subclass (the lone disagreement is 
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Branch et al. ( 2009 1 is the relative number of objects in each 
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Figure 10. The pEW of Si II A5972 versus the pEW of Si II 
A6355 of 89 SNe. Colours and shapes of data points are the same 
as in Figure [4] The solid lines are the boundariess between the 
four classes from [Branch e t al. (2009): CN = core normal, BL 
= broad line, CL = cool, SS = shallow silicon. The dotted lines 
are the median values of the two features' pEWs. The median 
uncertainty in both directions is shown in the upper-right corner. 



SN 1999ac, which is a borderline case between CN and SS). 
Classifications for the SNe shown in Figure [lO] can be found 
in Table [Al] in the "Branch Type" column. 

A slightly larger range of spectral ages than what was 



used by Branch et al. ( 2009 I is adopted here in order to in- 



crease the number of objects studied. The classifications and 
general trends seen when using only spectra within 3 d of 
maximum brightness are preserved when spectra within 5 d 
of maximum are used. Conversely, if a larger range of ages is 
used, the pEWs of the Si II A5972 and Si II A5972 features 
begin to evolve noticeably. The trends seen in Figure [To] are 
also visible at similar levels when plotting ApEW values 
(instead of pEW values), but by definition ApEW relies on 
a fit to the measurements as opposed to the measurements 
themselves. Finally, when the BSNIP dataset contains mul- 
tiple spectra of a given object within 5 d of maximum, only 
the spectrum nearest maximum brightness is used. This also 
has little to no effect on the classification of any SNe, nor 
on the trends seen in Figure [l0| 

One distinction between the two studies is that the 
BSNIP sample is lacking some of the most extreme members 
of the SS subclass. Most of these objects are the extremely 
peculiar SN 2000cx and SN 2002cx-like objects which are 
ignored in this study a priori. Other SNe in the bottom-left 
corner of Figure 2 of |Branch et al.| ( |2009| ) are Ia-91T, and 
even though there are a handful of Ia-91T SNe in the dataset 
presented here, the BSNIP sample does not contain spectra 
of any of these objects within 4 d of maximum. It does, 
however, include four Ia-99aa objects plotted in Figure [To] 
(SNe 1998es, 1999aa, 1999dq, 2001eh) which are thought to 
be intermediate objects between Ia-91T and la-norm (e.g., 
Garavini et al. 20041. They are much less clustered than 



class. Even though the total number of SNe is comparable in 
the two datasets, the BSNIP sample has significantly more 
CL and BL objects, mostly at the expense of CN objects (if 
one accounts for the fact that we are biased against some of 
the extreme SS SNe as mentioned above). 

One of the most striking trends seen in Figure [10] is 
that all of the data seem to form a nearly continuous, fairly 
well-correlated distribution with the most spectroscopically 
peculiar objects lying at the extreme edges of the distribu- 
tion. We have no Ia-91T objects in the figure, but they would 
fall in the bottom-left corner. From there, as the pEW of 
Si II A6355 increases at approximately constant Si II A5972 
pEW, we find the Ia-99aa, then la-norm, and finally HV 
SNe. Similarly, if we start at the median pEW values for 
the two features (where the dotted lines cross in Figure |l0| ), 
and increase Si II A5972 pEW at constant Si II A6355, we 
first find (mostly) la-norm and then the Ia-91bg objects. 
This amount of separation between SN la subclasses is not 
seen when comparing the pEWs of any other pair of features 
measured here. However, a somewhat weaker version of this 
separation is seen when comparing the pEW of Si II A6355 
to the pEW of both Ca II H&K and the S II "W." In ad- 
dition, Ia-91bg objects are found to distinguish themselves 
from the other subclasses mentioned here when comparing 
pEW values from a few pairs of spectral features. 

The subtypes used here (excluding HV) come from 
SNID, which cross-correlates the entire input spectrum with 
a library of SN spectra of various subtypes. However, it ap- 
pears that we can classify the vast majority of SNe la using 
only the pEWs of these two Si II features, as opposed to a 
spectrum covering a much wider wavelength range. The fact 
that there seems to be a continuous distribution in these 
pEW values with the "most" spectroscopically peculiar ob- 
jects on the edges of the distributiorj^] indicates that the 
often-used spectral classification scheme based on Ia-91bg, 
la-norm, and Ia-91T objects may only accurately represent 
the most extreme objects. In other words, an object spec- 
trally classified as la-norm might in reality have some ob- 
servables in common with la-norm and some in common 
with one of the peculiar subtypes. This classification scheme, 
while being somewhat qualitative in nature, is still useful, 
however, since any object spectroscopically classified as pe- 
culiar will be an outlier in the population of all SNe la, 
and it is often by studying the most extreme cases that one 
learns the most. A deeper, quantitative investigation into the 
photometric and host-galaxy properties of these spectroscop- 
ically determined subclasses will be undertaken in BSNIP III 
and other papers in this series. 



5.6 The Si II Ratio 

The so-called "Si II ratio,' 



K(Si II), was defined by |Nugent| 
et al. (119951) as the ratio of the depth of the Si II A5972 



they are in Branch et al. (20091, even though SN 1998es 



and SN 2001eh lie nearly on top of each other in Figure [10| 
The other major difference between this sample and that of 



5 In fact, most of the Ia-91bg and all of the Ia-99aa objects that 
appear in Figure [To] were used as SNID templates in BSNIP I. 
This indicates that they are indeed some of the "most" spec- 
troscopically peculiar objects in the BSNIP dataset, and that 
their spectra most closely resemble their subclass' namesakes 
(SN 1991bg and SN 1999aa, respectively). 
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Figure 11. The Si II ratio defined using spectral feature depth 
versus the Si II ratio defined using pEWs for 162 SNe. Colours 
and shapes of data points are the same as in Figure [4] The solid 
line is the linear least-squares fit and the dashed lines are the 
standard error of the fit. The Spearman rank coefficient is ~0.86. 



feature to the depth of the Si II A6355 feature. In our no- 
tation this is a(Si II A5972) / a(Si II A6355). They present 
a spectroscopic and photometric sequence of SNe la based 
on temperature differences which they attribute to differ- 
ences in the total amount of 56 Ni produced in the explo- 
sion. The Si II ratio has been previously found to corre- 
late with absolute _B-band magnitude (|Nugent et al.||1995[ ), 



Amis (|Benetti et al. 



2005l|Hachinger et al.|2006[ ), (B - V) Q 



( Altavilla et al.|2009 1 , and colour-corrected Hubble residual 
(Blondin et al. 112011 1. We will explore some of these rela- 



tionships in BSNIP III. 



While Nugent et al. (19951 originally defined the Si II 



ratio using the depths of spectral features, |Hachingcr et al~| 
( |2006| > define the Si II ratio using the pEWs of the Si II 
A5972 and Si II A6355 lines. With their redefined Si II ratio, 



Hachinger et al. ( 2006 1 observe a similar trend to what was 
found by |Nugent et al. ( |1995[ ) : brighter SNe la with broader 
light curves have smaller Si II ratios. In Figure [TT] we inves- 
tigate how consistent the Si II ratio is when defined using 
spectral feature depths (ordinate) versus pEWs (abscissa). 
The linear least-squares fit to all 162 points is shown as a 
solid line while the standard error of the fit is shown as the 
dashed lines. 

The two different ways to define the Si II ratio appear 
to be well correlated; the Spearman rank coefficient is '^0.86. 
According to the fit, the zero-point offset between the two is 
~0.06, though this should formally be zero since if the pEW 
is zero then the feature's depth is also zero. To measure the 
spectral feature depth one must determine the minimum of 
the spectral feature (either by smoothing the data or fitting 
a function to the flux, in order to avoid local, unphysical 
minima due to noise) and then define the feature's endpoints 
and a pseudo-continuum to compare to the flux at the mini- 
mum. However, to measure a pEW, one only needs to define 
the endpoints and a pseudo-continuum. Therefore, the pEW 
is a more robust and easier to measure parameter than the 
spectral feature depth, and since they are well correlated we 



Figure 12. The Si II ratio versus rest-frame age relative to maxi- 
mum brightness of 200 spectra of 154 SNe la. Colours and shapes 
of data points are the same as in Figure [4] 



will follow Hachinger et al. ( 2006 1 and define the Si II ratio 

as 



K (Si II) 



pEW (Si II A5972) 
pEW (Si II A6355) ' 



(3) 



The temporal evolution of the Si II ratio is shown 
in Figure |12| where we have once again colour-coded the 
data based on their HV or normal- velocity classification and 
shape-coded the data based on their SNID classification. 
There is a significant amount of scatter in 5R(Si II) and thus 
it is difficult to discern whether there is much temporal evo- 
lution for any of the subclasses of SNe la. The one point 
that can be fairly robustly made from inspecting Figure fl2] 
is that Ia-91bg objects consistently have the largest Si II ra- 
tios. This partially matches what has been previously seen 
in that fainter SNe la with narrow light curves (oftentimes 
spectroscopically resembling SN 1991bg) tend to have large 
jft(Si II) values (e.g.,|Nugent et al.|1995||Benetti et al.|2005 



|Hachinger et al.|2006 |. However, these previous results also 
show that more luminous objects with broader light curves 
have the smallest Si II ratios, which is not readily appar- 
ent in the current study. There are only a few Ia-91T/99aa 
SNe plotted in Figure |12| but they all have average values 
of K(Si II). 

Figure [12] also suggests that HV objects may have, on 
average, lower Si II ratios than SNe with normal velocities, 
although the large scatter renders this conclusion tentative 
at best. When colour-coding the data by "Benetti type," 
as expected the FAINT objects have the highest values of 
5ft(Si II) (since they are defined as SNe having the narrow- 
est light curves) . Somewhat different than previous studies is 
that LVG and HVG objects all have approximately the same 
values of the Si II ratio. Furthermore, while we only have a 
handful of objects with a measured 5R(Si II) at t < —5 d, a 
few are LVG and a few are HVG, and they all have approx- 
imately average values of the Si II ratio. This is markedly 
different than reported by |Benetti et al.| ( |2005[ l, whose data 
show that HVG objects have larger than average 5R(Si II) 
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Figure 13. The pEW of the Si II A5972 feature versus the Si II 
ratio for 89 SNe la. Colours and shapes of data points are the same 
as m Figure[I] The various lines are fits to all of the data (solid), 
only la-norm objects (dotted), only HV objects (dot-dashed), 
only Ia-91bg objects (long-dashed), and only Ia-91T/99aa objects 
(short-dashed). All of the fits have Spearman rank coefficients of 
> 0.89 except the Ia-91bg objects (which have a coefficient of 
~0.79). 



values at these earliest epochs (though they do decrease to 
more typical values at later times). 

The pEW of the Si II A5972 feature and the Si II ra- 
tio have both been found to correlate well with Amis (e-g-, 
Nugent et al.|1995 Hachinger et al.||2006 |. Hachinger et al. 
( |2006| > point out (from their Fig. 10) that each of these 
correlations implies the other's existence based on the rela- 
tionship between the near-maximum values of pEW of Si II 
A5972 and 3?(Si II). In Figure p] the pEW of the Si II A5972 
feature versus the Si II ratio is plotted for objects within 
5 d of maximum, confirming this correlation. The plot con- 
tains 89 SNe la. If the current dataset had multiple spectra 
of a given object for which 5R(Si II) and the pEW of the 
Si II A5972 feature was calculated, the observation closest 
to maximum brightness was usedQ 

In Figure 13 the current study's definition of 5R(Si II) 
is used (i.e 



the ratio of pEWs) as opposed to the origi- 
nal definition (i.e., the ratio of depths; Nug ent et al.|[T995| > . 
Figure 10 of Hachinger et al. ( 2006 1 uses the original defi- 



nition on the abscissa. When the original definition is used 
with the BSNIP data, results qualitatively similar to those 

~ (12006b interpo- 



of Figure 13 



Hachinger et al. 



late/extrapolate their pEW measurements to t = d, while 
spectra within 5 d of maximum brightness (as in Section [575| 
are used in Figure [l3| since the pEWs of the Si II A5972 and 
Si II A5972 features do not change much during these epochs. 
If the few objects with spectra earlier than 5 d before max- 
imum are included, they fall along the linear correlation. 
However, if objects with spectra later than 5 d after maxi- 
mum are added, they mainly fall below the linear correlation. 

6 This practice of using only the spectrum nearest to maximum 
brightness when multiple spectra of a given object are analysed 
is followed in the rest of this work unless otherwise noted. 



This is due to the fact that the Si II A5972 pEW values re- 
main relatively constant for t < 5 d but begin to increase at 
later times, as seen in Figure [7] (bottom left). 

As mentioned above, Ia-91bg objects have large Si II 
ratios and large Si II A5972 pEWs, and are thus found in 
the upper right of the linear correlation in Figure [13] Con- 
versely, it was shown that Ia-91T/99aa SNe have low values 
of Si II A6355 pEW and low-to-average values of Si II A5972 
pEW, so one may expect them to lie at the bottom left of 
the correlation, and in fact it appears that these objects lie 
somewhat below the main correlation. The HV SNe perhaps 
make up the upper part of the main correlation, with the 
normal- velocity objects making up the lower part, though 
these two groups are fairly well mixed. Similarly, if we code 
the points in Figure [13] by velocity gradient, we see that the 
HVG objects usually lie slightly above the LVG objects (as 
. Hachinger et al.|2006 1 , but with a significant amount 



of overlap. The FAINT objects in the BSNIP sample are 
found at the upper right of the correlation, as seen previ- 
ously (Hachinger ct al. 2006), and reinforce the connection 
between Ia-91bg SNe and FAINT SNe. 

To investigate whether the various subclasses have dis- 
tinct correlations, linear functions were fit to all of the data 
in Figure [13] as well as each individual subclass. All of the 
fits have Spearman rank coefficients of > 0.89 except the Ia- 
91bg objects (which have a coefficient of ~0.79), implying 
that each subclass, and the data as a whole, are well fit by 
a linear function (though Ia-91bg SNe have a fair amount of 
scatter). The linear fit to all of the data and the linear fits 
to the la-norm, HV, and Ia-91bg objects are all consistent 
with each other at the 2a level. As mentioned above, the Ia- 
91T/99aa SNe fall below the main correlation and appear to 
lie on their own relationship. However, since there are only 
four points in the fit, it is difficult to say if this difference 
is statistically significant. Previous work has indicated that 
both the pEW of the Si II A59 72 feature and the Si II ratio 
are luminosity indicators (e.g., Hachinger et al.| [2006). The 
scatter in the relationship between these two parameters 
may be related to differences in luminosity or light-curve 
shape (i.e., Ia-norm versus Ia-91bg versus Ia-91T) or dif- 
ferences in velocity (i.e., HV versus normal velocity), which 
might also be related to differences in colour ( Foley & Kasen 



20111. Adding photometric observables such as these into 



this analysis will be done in BSNIP III. 

The pEW of the Si II A4000 feature has also been used 
as a luminosity indicator |Arsenijevic et al.||2008| |Walker| 



et al.|2011||Blondin et al.|2011||Chotard et al.|201l| . There 



fore, one might expect it to also correlate well with the Si II 
ratio. Figure [14] shows the pEW of the Si II A4000 feature 
versus 5R(Si II) for objects within 5 d of maximum bright- 
ness. The plot contains 66 SNe la. 

Figure [14] shows no obvious correlation between pEWs 



of the Si II A4000 feature and 5R(Si II). As in Figure [13 
Ia-91bg are found to have the highest Si II ratio and are 
thus toward the right of the figure, while the Ia-91T/99aa 
SNe have small pEWs in general and are thus found to- 
ward the bottom of the plot. The HV objects have higher- 
than-average Si II A4000 pEW values with relatively low 
Si II ratios, but they overlap quite a bit with the normal- 
velocity objects. The Ia-norm SNe show the most consistent 
relationship between these two parameters, although they 
have a large scatter and a Spearman rank coefficient of only 
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Figure 14. The pEW of the Si II A4000 feature versus the Si II 
ratio for 66 SNe la. Colours and shapes of data points are the 
same as in Figure [4] 



Figure 15. The ratio of the pEW of S II "W" to that of Si II 
A6355 versus SiS ratio for 110 SNe la. Colours and shapes of data 
points are the same as in Figure [4] 



~0.51. The fact that these two spectroscopic observables are 
both used as luminosity indicators (e.g., Blondin et al.|20TT I 
and yet they are seemingly uncorrelated is curious. Further 
investigation into this discrepancy (with the addition of lu- 
minosity information for some of the SNe presented here) 
will take place in BSNIP III. 

5.7 Other pEW and Flux Ratios 

5.7.1 The Call Ratio 



Similar to the Si II ratio, Nugent et al. ( 1995 I defined the 



Ca II ratio as the ratio of the flux at the red edge of the 
Ca II H&K feature to the flux at the blue edge of that fea- 
ture. In our notation this is defined as 

F r (Ca II H&K) 



5ft (Ca II) 



(4) 



F b (Ca II H&K) ' 
They found that this parameter scaled with absolute B-band 
magnitude in the same way as 5ft(Si II). The Ca II ratio 
has been investigated further in more recent studies, and 
while there is not much temporal evolution of the parameter 
within 5 d of maximum brightness ( Bongard et al.|2006 l, it 
is unclear whether this spectral parameter is reasonably well 



correlated with luminosity ( Blondin et al.|2011 l. 

Since the Si II ratio and the Ca II ratio have both been 



shown to correlate with absolute B-band magnitude (Nu- 
gent et al.|1995| |Bongard et al.|2006| ) , one might expect the 
two ratios to correlate with each other. However, we find no 
significant correlation. 

5.7.2 The "SiS" Ratio 



Somewhat analogous to the Ca II ratio, |Bongard et al.| 
( 2006 1 defined the "SiS ratio" as the ratio of the flux at 



the red edge of the S II "W" feature to the flux at the red 
edge of the Si II A6355 feature. This is defined using our 
notation as 



5ft (SiS) 



F r (S II "W") 
F r (Si II A6355) ' 



(5) 



This parameter was shown to scale with absolute B-band 



magnitude in the same way as 5ft(Ca II). Even though Bon- 



gard et al. ( 2006 1 make a time-dependent correction to yield 



this correlation, they state that the existence of a universal 
time correction for this ratio is doubtful. In fact, very little 
temporal evolution of 5ft(SiS) near maximum brightness is 
seen in the BSNIP data. The SiS ratio has been measured for 
more objects in recent work, but it does not appear to cor- 



relate as well with luminosity as initially claimed (Blondin 
et al.|2011 1 



Since the SiS ratio and the Si II ratio have both been 



shown to correlate with absolute B-band magnitude (Bon 
|gard et al-lpOo) , one might expect the two ratios to cor- 
relate with each other. Effectively no evidence for such a 
correlation is found in the current dataset. 

One may expect the SiS ratio to be related to the pEW 
of the S II "W" feature or the Si II A6355 feature (or per- 
haps the ratio of the two). We find no evidence for such a 
correlation with the pEW of the S II "W" alone, but a pos- 
sible relationship exists between 5ft(SiS) and the pEW of the 
Si II A6355 feature, as well as the ratio of the pEW of the 
S II "W" to that of Si II A6355. Figure [15] presents 110 SNe 
within 5 d of maximum brightness for which we measure 
5ft(SiS). Increasing or decreasing the age range investigated 
does not change the basic trends seen in Figure [15| 

It is possible that the SiS ratio is linearly (or quadrati- 
cally) related to the ratio of pEWs, but there is a significant 
amount of scatter in Figure |15| Interestingly, Ia-91bg and 
Ia-99aa objects once again tend to lie on the outskirts of 
the distribution, but there is much overlap with the HV 
and la-norm objects (which are nearly indistinguishable in 
this parameter space). There are also a few extreme out- 
liers in Figure [l5| including the lone Ia-91T object plotted. 
When comparing 5R(SiS) to the pEW of Si II A6355 alone, 
the same separations (or lack thereof) between the different 
subclasses once again appear, as do the same outliers. The 
only difference is that the two values are (possibly) linearly 
anticorrelated. 
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Figure 16. The Si II ratio versus the SSi ratio for 85 SNe la. 
Colours and shapes of data points are the same as in Figure [4] 
The solid line is a linear fit to all of the data and the dashed line 
is a linear fit to only the la-norm and HV objects. The Spearman 
rank coefficient for both fits is approximately —0.89. 



5.7.3 The "SSi" Ratio 

Another parameter that has been proposed as a spectro- 
scopic luminosity indicator is the ratio of the pEW of the 
S II "W" to that of the Si II A5972 feature (Machinger et al. 



2006). They found that this ratio, dubbed 5R(S,Si), is lin- 



early anticorrelated with Amis (which is opposite the rela- 
tionship between 5ft(Si II) and Amis). Following Hachinger 
et al. (20061, we define the "SSi ratio" as 



SR(S,Si) 



pEW (S II "W" 



(6) 



pEW (Si II A5972) ' 

There is strong evidence for an anticorrelation between the 
Si II ratio and the SSi ratio in Figure [To] which presents the 
85 SNe (within 5 d of maximum) for which both 3?(S,Si) and 
5R(Si II) are measured. 

Interestingly, the overall correlation (and differences be- 
tween different subclasses) is unchanged when we include 
data from all epochs studied in this work. It is also un- 
changed when a narrower age range is used. The overall 



anticorrelation seen in Figure 16 suggests that 5R(S,Si) may 
indeed be as good a luminosity indicator as K(Si II), though 
recent work suggests otherwise ( Blondin et al.|2011 l. While 
the trend for all objects shown in Figure |16| may not be lin- 
ear, they do appear well correlated, with a Spearman rank 
coefficient of about —0.89. A linear fit to all 85 SNe in Fig- 
ure [16] is shown by the solid line. A linear fit to only the 
la-norm and HV objects (shown by the dashed line in Fig- 
ure 161 has a nearly identical Spearman rank coefficient to 



that of the entire sample. The tantalising possibility that 
5R(S,Si) is yet another spectral luminosity indicator will be 
explored further in BSNIP III. 

5.7.4 The "SiFe" Ratio 



Similar to the SSi ratio, Hachinger et al. (2006) also found 



that the ratio of the pEW of the Si II A5972 feature to that of 
the Fe II complex can be used as a spectroscopic luminosity 



Figure 17. The Si II ratio versus the SiFe ratio for 72 SNe la. 
Colours and shapes of data points are the same as in Figure [4] 
The solid line is a linear fit to all of the data and the Spearman 
rank coefficient is ~0.86. 



indicator since it scaled linearly with Amis. They referred 
to this as the "SiFe ratio," defined as 



5R (Si,Fe) 



pEW (Si II A5972) 
pEW (Fe II) ' 



(7) 



There is also strong evidence for a correlation between 
the Si II ratio and the SiFe ratio in Figure |17| where the 
72 SNe (within 5 d of maximum) are shown for which both 
SR(Si,Fe) and 3?(Si II) are measured. 

The scatter in the relationship between the Si II ra- 
tio and the SiFe ratio increases when including data from 
all epochs studied here. The overall correlation seen in Fig- 
ure 17 suggests that 5R(Si,Fe) may also be an accurate lu- 



minosity indicator, though Blondin et al. (2011 1 report that 



it does not significantly improve the accuracy of their lumi- 
nosity determinations. The data are well fit by a linear func- 
tion and are highly correlated (Spearman rank coefficient of 



4). 86); a linear fit to all 72 SNe in Figure 16 is shown by the 



solid line. The Ia-91bg objects distinguish themselves quite 
prominently except the lone point that appears to be closer 
to the Ia-99aa objects, which are also found at an outer 
edge of the main relationship. The HV SNe tend to be at 
the bottom of the relationship, with a handful of significant 
outliers. The use of 5R(Si,Fe) as a luminosity indicator will 
be investigated further in BSNIP III. 



5.8 



Comparing Expansion Velocities to 
Pseudo-Equivalent Widths 



In Section [5. 4. 7| the possibility of a correlation between ex- 
pansion velocity and pEW was mentioned, at least for the 
Ca II H&K, Mg II, Fe II, and Si II A6355 features. The pEW 
for these four features is plotted versus the expansion ve- 
locity of Si II A6355 in Figure [18] They contain all objects 
within 5 d of maximum brightness. In all cases, decreasing 
the age range investigated does not change the basic trends 
seen in the figure, though increasing the age range increases 
the overlap among the various subclasses. 
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Figure 18. The pEW of various spectral features versus the velocity of Si II A6355. (top left) 92 SNe la for Ca II H&K, (top right) 99 
SNe la for the Mg II complex, (bottom left) 99 SNe for the Fe II complex, and (bottom right) 120 SNe for Si II A6355. Colours and shapes 
of data points are the same as in Figure[4] The median uncertainty in both directions is shown in the lower-right corner of each panel. In 
the bottom-right panel the filled blue square is the HV object SN 2006X which shows evidence for interaction with circumstellar material 
(Patat ct al. 2007f and an aspherical explosion ( |Patat et al.|2009| , and the filled black squares are the HV objects SNe 2002bf, 2002bo, 
and 2004dt, all of which show evidence for an aspherical explosion via large Si II polarizations ||Wang et al.|2006). 



By construction, HV objects are always found toward 
the right-hand side of plots of pEW versus Si II A6355 ve- 
locity. But, as hinted in Section [5. 4. 7| they also have higher- 
than-average pEW values for most of the spectral features 
investigated. The features showing the greatest differenti- 
ation between la-norm and HV objects in this parameter 
space are the four mentioned above and plotted in Figure [l8l 

In each of these four Figures there is significant overlap 
among spectroscopically peculiar objects and la-norm, while 
the HV SNe tend to be fairly well separated. However, there 
is a large scatter, and thus it is difficult to say whether there 
is any real correlation between pEW and Si II A6355 velocity 
for the la-norm and the HV objects. The Spearman rank 
coefficients for the panels in Figure [l8| are all about 0.21- 
0.35. Using only the la-norm and HV objects, the Spearman 
rank coefficients improve only moderately for these spectral 
features (0.28-0.33). 

The bottom-right panel of Figure [18] shows the largest 



difference between the various subclasses of SNe la, although 
there is still significant overlap (mostly between la-norm and 
Ia-91bg objects). There is a more significant correlation be- 
tween the pEW of Si II A6355 and the expansion velocity 
derived from that same feature, though it is still not too 
significant (Spearman rank coefficient of 0.43 for all objects 
in the bottom-right panel of Figure [Ps] and 0.51 for la-norm 
and HV SNe only). The basic trend seen in this figure has 



also been noted previously by Wang et al. ( 2009 1 with some 



of the same spectra used in this study. Their plot shows only 
SNe la within 3 d of maximum while Figur e |18| p resents ob- 
jects within 5 d of maximum. Wang et al. (2009) point out 



that HV objects typically have pEW values greater than 
~100 A, which is also observed here. 

The observation of a strong Si II A6355 feature at high 
velocities may indicate an enhancement in the abundance 
or density in the outermost layers of HV objects. This could 
be caused by an extended burning front, different degrees 
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of mixing in the ejecta, or interaction with circumstellar 
material (CSM; |Benetti et al.|[2004| |Wang et al.||2009[ ). In 

fact, 



two SNe la with some of the best direct observational 
evidence for CSM interaction, SN 2006X ( |Patat et al.|2007[ 
the filled blue square in the bottom-right panel of Figure |18p 
and SN 20071e ( jSimon et aL||2009[ ), are both in the current 
sample and are found to be HV objects. 



Wang et al. ( 2009[) also point out that the distribution 



of pEWs and velocities between la-norm and HV objects 
is continuous, and therefore it is likely that the proposed 
enhancement occurs at a wide range of values (with the ma- 
jority of SNe la showing relatively low levels of this enhance- 
ment). They suggest that this could be a result of viewing 
angle effects if HV objects are associated with asymmetric 
structures. There is both observational (Leonard et al.|2005 



Wang et ai]|2006| > and theoretical ( |Maeda et al.||20ld[ ) evT- 



dence that HV SNe may be associated with certain viewing 
angles of aspherical explosions. 

Leonard et all (120051) , IWang et all 12007L and 



Patat et all (I2009I) discuss spectropolarimetric data on 



SNe 1997bp, 2002bf, 2002bo, 2004dt, 2004ef, and 2006X, all 
of which show large Si II polarizations, implying an asym- 
metric explosion. These objects are found to be HV SNe in 
the BSNIP data; SNe 2002bf, 2002bo, and 2004dt are plotted 
as filled black squares in the bottom-right panel of Figure fl8] 
and SN 2006X is plotted as a filled blue square. Building on 



the work of Leonard et al. (20051, Maund et al. (20101 found 



a correlation between the level of Si II polarization and v. Of 
the six aforementioned HV SNe with large polarizations, the 
two for which we measure a velocity gradient are both HVG 
objects. Si II A6355 velocity, v, and level of polarization may 
all also be correlated with velocity offsets of nebular lines in 
late-time spectra ( Maeda et al.|20To I . Comparisons of spec- 
tra at later epochs to models such as those presented by 



Maeda et al. (20101 will be made in future BSNIP studies. 



6 CONCLUSIONS 

This is the second paper in the BSNIP series, and it presents 
spectral feature measurements of 432 low-redshift (z < 0.1) 
optical spectra within 20 d of maximum brightness of 261 
SNe la. These data all come from BSNIP I (Silverman et al, 
submitted) and were obtained from 1989 through the end of 
2008. We have outlined in detail our automated and robust 
procedure for spectral feature definition and measurement 
which has expanded upon the work of previous studies. Us- 
ing this algorithm we attempted to measure expansion ve- 
locities, pEWs, spectral feature depths, and fluxes at the 
centre and endpoints of each of nine major spectral feature 
complexes. The raw numbers measured from the data are 
also presented. A sanity check of the consistency of our mea- 
surements was performed using the BSNIP data in addition 
to the spectra presented by Mathcson et al.| ( |2008[ ). Even 
though the current study utilises spectral data almost ex- 
clusively, future BSNIP papers will incorporate photometric 
data and host-galaxy properties into the measurements dis- 
cussed in this work. 



6.1 Summary of Spectral Feature Measurements 

The temporal evolution of the expansion velocity of each of 
the spectral features (except for the Mg II and Fe II com- 
plexes) was explored. The observed differences in velocities 
support the layered structure of SN la ejecta, with Ca II 
found in the outermost (i.e., fastest expanding) layers and 
O I, Si II, and S II found in the inner (i.e., slower expand- 
ing) layers. Many of the basic trends seen previously are 
confirmed by the current analysis, including evidence of a 



HV population of spectroscopically normal SNe (Wang et al 



2009). A sharp cutoff between normal- velocity and HV ob- 



jects is not well motivated by the data analysed here, though 
there are some differences between spectroscopically normal 
SNe la with the lowest and highest velocities near maximum 
brightness (which is discussed in more detail below). 

Following |Benetti et al. (20051, velocity gradients are 



calculated and then used (along with light-curve shape, for 
those that have that information) in order to classify the 
SNe. The BSNIP dataset is not the best suited for this kind 
of study due to the relatively low number of spectra per 
object. However, v is still measured for 61 objects. 

The interpolated/extrapolated velocity at maximum 
brightness (vo) and at 10 d past maximum (i>io) was cal- 
culated for each SN with a measured i>. Effectively no corre- 
lation is found between vio and velocity-gradient classifica- 
tion, and only a weak correlation is seen between vo and ve- 
locity gradient. In previous work HV and HVG objects have 
been used almost interchangeably, as have normal-velocity 
and LVG objects (e.g., Hachinger et al.|2006 Pignata et al 



2008 



Wang et al.||2009 l. However, the data presented here 



cast serious doubt on these associations. 

The temporal evolution of the pEW of each spectral fea- 



ture was also examined. As in Nordin et al. (2011b I, ApEW 



was calculated by subtracting a fit to the data from the 
measurements themselves. Unlike previous work, only lin- 
ear and quadratic functions were used to fit the temporal 
evolution. ApEW is supposed to effectively remove the age 



dependence on pEW measurements (Nordin et al. 2011b I, 



though the scatter in the pEWs of most features is large 
enough in the BSNIP data that there is still some level 
of temporal evolution in ApEW. Furthermore, the ApEW 
values rely on defining a fit (either linear or quadratic) to 
the measurements, adding another assumption to the anal- 
ysis. No correlation between velocity gradient and pEW or 
ApEW is seen in the data presented here, even though there 
has been evidence presented previously for a relationship be- 



tween ApEW of the Si II A4000 feature and v (Nordin et al 



2011a) 



We use pEW values of Si II A6355 and Si II A5972 to 
classify the BSNIP data according to the four groups de- 



fined by Branch et al. ( 2006 1. The boundaries between these 



groups appear to be arbitrary in the context of the spectral 
measurements presented here. That being said, spectroscop- 
ically peculiar SNe la and HV objects occupy the outermost 
edges of the nearly continuous parameter space filled in by 
the la-norm objects. Thus, the "SNID types" seem to be the 
most extreme versions of the three non-CN "Branch types." 
Furthermore, both of these classification schemes are quasi- 
arbitrary cuts in what appears to be a continuous distribu- 
tion of pEW values. 

It was shown that the Si II ratio can be equivalently 



BSNIP II: Initial Spectral Analysis 27 



defined as a ratio of spectral feature depths or pEW values 
and the latter is used as the definition for the present study. 
This quantity has been shown to correlate with a variety of 



photometric observables (e.g., Nugent et al. 1995 



et al. 2006 Altavilla et al. 2009 Blondin et al 



Hachinger 



2011) and 



these proposed relationships will be investigated further in 
BSNIP III. We do confirm the observation that Ia-91bg ob- 
jects (often underluminous) have the largest 5ft(Si II) values. 
On the other hand, the Ia-91T/99aa objects (often overlu- 
minous) seem to have average 5ft(Si II) values, which differs 
from what has been previously seen (e.g., Hachinger et al. 
20061). 



Benetti et al. ( 2005 I presented evidence of HVG objects 
having the largest measured Si II ratios at t < —5 d, though 
the few objects for which we measure 5R(Si II) at these early 
epochs all have quite typical Si II ratios and represent the 
LVG and HVG classes. The pEW of the Si II A5972 feature 
and 5R(Si II) have both been found to correlate well with 
Amis (e.g., Nugent et al. 1995; Hachi nger et al.|2006 l, and 
the BSNIP data show strong evidence of a linear correlation 
between these spectral parameters. While the HV and Ia- 
norm objects show quite a bit of overlap in this parameter 
space, the Ia-91bg objects lie at the uppermost end of the 
correlation and the Ia-91T/99aa objects are well below the 
main trend (and perhaps form there own, distinct relation- 
ship). 

Four other pEW and flux ratios were also calculated 
and discussed. Two of these, the Ca II ratio and SiS ra- 
tio, have been previously seen to be correlated with B-band 



magnitude in the same way as the Si II ratio ( Nugent et al. 
1995 Bongard et al.|[2006 respectively). However, no rela- 
tionship between 3i(Si II) and 5ft(Ca II) is seen here, and only 
a weak relationship is observed between 5R(Si II) and 3?(SiS). 
The various "SNID types" and 'Wang types" all have sim- 
ilar ranges of Ca II and SiS ratios. The other two ratios 
investigated, the SSi ratio and the SiFe ratio, were shown 



by Hachinger et al. ( 2006 1 to be anticorrelated and corre- 



lated with Amis, respectively. There is strong evidence for 
both of these trends in the BSNIP data. Most of the "SNID 
types" and "Wang types" have a high degree of overlap in 
the plots of Si II ratio versus SSi ratio and SiFe ratio, and 
the Ia-91bg objects occupy the extreme upper end of both 
of these correlations. 

Finally, possible correlations between the expansion ve- 
locity of Si II A6355 and pEWs of Ca II H&K, Mg II, Fe II, 
and Si II A6355 were explored. Spectroscopically peculiar 
objects mostly overlap with la-norm objects in these plots, 
while the HV SNe show both higher expansion velocities and 
pEW values. CSM interaction and asymmetric explosions 
have been proposed to explain the existence of SNe la with 
high velocities and large pEWs ( |Benetti et" al. 2005 |W"ang| 
et al.|2006 respectively). Interestingly, four of the points at 
the uppermost end of the Si II A6355 velocity-pEW relation- 
ship show evidence for CSM interaction (SN 2006X) or an 
aspherical explosion (SNe 2002bf, 2002bo, and 2004dt). 

6.2 How Should One Spectroscopically Classify 
SNe la? 

Throughout this paper, various SN la spectral classification 
schemes have been used, but which one is "best" or the most 
illustrative? The "Benetti types" based on velocity gradient 



measurements are not only difficult to determine, requiring 
velocity measurements of multiple spectra near maximum 
brightness, but the various subclasses do not show many ob- 
vious observable differences (aside from v and Amis values, 
which is how the groups are defined). As mentioned above, 
the "SNID types" are merely the most extreme objects of 
the non-CN "Branch types," though the delineation between 
the different subclasses is fuzzy at best. Finally, the distinc- 
tion between HV and normal- velocity objects is somewhat 
arbitrary as well. 

Despite the ambiguity in subclass determination, we 
prefer to use the "SNID types" (la-norm, Ia-91bg, Ia-91T, 
and Ia-99aa) as defined in BSNIP I. We also favor split- 
ting the la-norm objects into HV and normal-velocity ob- 
jects as defined in Section [5. 1| It should be made clear that 
these distinctions may merely be a naming convention used 
for the most extreme SNe la (according to some measured 
observable), and in reality all of these objects represent a 
continuous distribution of observables. 

One can ask, "How 91bg-like is an object?" if measure- 
ments of some spectral parameters are similar to Ia-91bg 
objects and some parameters are more similar to la-norm. 
We suggest that this is an ill-posed question and that one 
should be less concerned with naming conventions and more 
concerned with how the various spectral measurements cor- 
relate with luminosity (or other observables). One might re- 
spond to this with the question, "Why keep the naming con- 
ventions at all?" The spectroscopically peculiar SNID types 
and the HV objects are the extreme members of the SN la 
population and it is often by studying the most extreme ob- 
jects that one understands the more typical specimen. This 
is why we prefer the classification scheme outlined above. 

6.3 Can a Theoretical Model of SNe la Explain 
Their Spectra? 

Many of the spectral measurements investigated in this 
study show a continuous (or nearly so) range of values. This 
is perhaps evidence that all of the SN la subclasses consid- 
ered here can be described by one, self-consistent physical 
model that only varies in one or two intrinsic parameters. 
It was suggested well over a decade ago that the Ia-91bg, 
la-norm, and Ia-91T/99aa objects can be naturally linked 
through a continuous increase in temperature, presumably 
by a smoothly increasing amount of 56 Ni produced in the 
explosion (e.g., N ugent et al.|[l995[ ). More recently, it has 
been suggested that continuous variations in viewing angle 



can give rise to HV and normal- velocity objects (Kasen & 
|Plewa|[2007l |Kasen et aL][2009l |Maeda et aT]|2010[ ). Any" 
model that hopes to reflect reality must match the observed 
spectral properties near maximum brightness (as presented 
here), as well as the photometric observations of SNe la (as 
shown by Ganeshalingam et al. 2010[ and discussed further 
in BSNIP III). 

Perhaps by varying both the amount of 56 Ni produced 
and the viewing angle, one can self-consistently explain all 
of the subclasses of SNe la (and their measured spectral pa- 
rameters). Though tantalizing, a deeper investigation of this 
possibility is beyond the scope of this paper. Furthermore, if 
one (or more) of the subclasses of SNe la are in reality funda- 
mentally different, then discussing only a single theoretical 
model does not make sense. At least two distinct progenitor 
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channels (possibly resulting in different observed spectra) 
have been suggested for SNe la: the single degenerate (e.g., 



Whelan fe Iben|1973[) and t he double degenerate (e.g., Iben 
fe Tutukov||1984||Webblnkl[T984| . Similarly, different types 



of SNe fa may be found in different host galaxies (perhaps 
related to host-galaxy type, mass, metallicity, or stellar age; 



e.g., [Gallagher et al-pOOH] |Hicken et al.|[2009l |Howell et aL] 



20091, which may also point to multiple models being re- 



quired to explain all of the observed spectra. 



6.4 What About Future SN Surveys? 

Surveys which are much larger than BSNfP and tuned to 
discovering the bulk of their SNe la at higher redshifts than 



BSNIP are already underway (e.g., Rau et al. 2009 Law 



ct al. 2009; |Kaiser et al.|2 002), with many more planned for 



the future (e.g., LSST, WFIRST). With such large numbers 
of objects being discovered, it will be extremely difficult to 
obtain high-quality light curves and/or multiple spectra of 
most SNe. fnstead, one will need to rely on only a handful 
of photometric points and perhaps one, often relatively low 
S/N, spectrum. Thus, spectral indicators of luminosity that 
can be measured in one low-quality spectrum will be of the 
utmost importance. 

For higher-z surveys in particular, useful spectral fea- 
tures that are found toward the blue end of the optical range 
will be even more critical. For z > 0.6, the red wing of 
the typical, near-maximum Si II A6355 feature becomes red- 
shifted beyond ~1 fim. The Ca II ratio is promising in this 
regard, though no relation between it and the Si II ratio was 
seen in the BSNIP data. The SiFe ratio does correlate well 
the Si II ratio, and the Fe II complex is easier to measure 
and bluer than the Si II A5972 feature, but it still requires a 
measurement of the pEW of Si II A6355. 

The accuracy as luminosity indicators of these and other 
spectral observables discussed herein will be explored in the 
next paper in this series, BSNIP III. It will utilise the mea- 
sured values and relationships described in this work, exam- 
ining the correlations between spectroscopic observables and 
photometric properties, with an eye toward using a single 
spectrum to determine the luminosity (and possibly other 
intrinsic characteristics) of a SN la. 
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APPENDIX A: SUMMARY OF SPECTRAL 
DATASET 

Table [AT] present s each SN la which had a measured pseudo- 
continuum for at least one spectral feature. Also listed is the 
(rest-frame) spectral age for each observation that was fit, 
and spectral classifications based on various classification 
schemes. 
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Table Al: Summary of Spectral Dataset 
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sn lyyygn 


a i fi 1 r on 1 cr ot 

4. 12, 15. 89, 15.97 


I a- norm 


IAIN 1 


CL 


N 


oat onnnL l 

biN 2000bk 


14.84 


I a- norm 








biN 200Ucn 


14.25 


I a- norm 








oat ooon„~-. 

SN 2000cp 


r» no 1 1 TCi 

2.92,11.70 


I a- norm 






N 


OAT OOOO 

SN ZUUUcu 


9.64 


I a- norm 






N 


oat ooon 

SN ZUOOcw 


4.81 


I a- norm 




BL 


N 


oat nnnoJ~ 

SN 2000dg 


c no a cr 

— 5.09,4.66 


I a- norm 




bb 


N 


oat innnJi 

biN 2000dk 


i on 1 n o a 

1.00,10.84 


I a- norm 


T7 A TATT 1 
IAIN 1 


CL 


N 


cm onnn^vn 
biN zUUUdm 


1 KQ fi 1Q 
— 1.00,0.10 


I a- norm 


tt"\ jr^ 
nVb 


u i 

L>Li 


IN 


OAT OOOOJ 

SN 2000dn 


O O A A f O O 

— 0.94,16.38 


I a- norm 


LVG 


BL 


N 


SN 2000dr 


6.78 


I a- norm 






N 


SN 2000dx 


— 9.26 


I a- norm 








SN 2000ey 


7.90 


I a- norm 






HV 


SN 2000fa 


-8.25,6.86 


I a- norm 






N 


SN 2001E 


15.01 


I a- norm 








SN 2001G 


11.57 


I a- norm 








SN 2001N 


13.05 


I a- norm 








SN 2001V 


15.86 


I a- norm 








SN 2001ay 


6.79 


I a- norm 






HV 


SN 2001az 


-3.24 


I a- norm 






N 


SN 2001bf 


1.22 


I a- norm 






N 


SN 200 lbg 


13.70,18.91 


I a- norm 


HVG* 






SN 200 lbr 


3.47,3.48 


I a- norm 




BL 


HV 
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Tabic Al — Continued 



SN Name 


Phase 3, 


SNID 


Bcnctti 


Branch 


Wang 






(Sub)Tvne b 


TvDe c 


TvDe d 


TvDe° 


SN 2001 hn 

Oil £t\J\J _L L/LI 


0.51 


±cX 11LJ1 111 








SN 2001 rn 


1 39 1 s oo 


± Cl llLJl 111 






N 


SN 2001 da 

Oil £i\J\J 1 Lid 


— 1.12.9.72 


±&< 11LJ1 111 


HVG 


BL 


N 


sn 2001 di 

Oil ^KJ\J ± 


13.84 


Tn-Tinvtn 

± CL LLVJL 111 








SN 2001 dt 

Oil iiL/\JH-ll/ 


13.60 


T a _n riTin 

± cX llLJl 111 








SN 2fim dw 


11.06 


Tn -Tinvm 

± cX llLJl 111 








SN 2001 ph 

i iiUVJlt_ll 


— 5 63 —4 53 3 26 


J. £X <X<Jj 




SS 




SN 2001 pn 

oi i £d\j\.J iv_.ii 


1 0Q 1 4 72 

1VJ . UCy , Irt . 1 


T a _ti nrm 
±cx iilji ill 


HVG 






SN 2nnien 

Oil Zj UVJ 1 "-'^J 


2 83 4 99 5 97 7 85 


Tn-Tinvtn 

±cX 11LJ1 111 


HVG 


CL 


N 


SN 2001 px 

Oil _-U\^l - A 


— 1.82 


ACL rj _L Utl 








SN 2fim fp 


— 0.99 


T a _ti Arm 

±CX 11LJ1 111 




SS 


N 


SN 2fini fh 

L)l "I __iUVJ 1 111 


5 93 7 84 


± CX11LJ1 111 






N 


sn 20020 

Oil __■ UVJZj \_T 


19.31 


T a _n riTin 
±cx iilji in 








SN 9nf)9aw 

Oil __rl_JVJ__idW 


9 i n 1 84 

— . J. \J , A U . Orfc 


T £» "Tll^VTYl 
id 11U1 111 






N 


SN 2fin2hf 


2 97 6 90 


±cX 11LJ1 111 




BL 


HV 


SN 9nf)9hn 

Oli -jUu-jUU 


— IT QA _1 OS 1 ^ QQ 


Ta Tij^T"Tn 

Id 11U1 111 


HVG 




HV 


SN 2002bz 


4.92 


Is,- norm 




CL 


N 


SN 2f)n2rd 


1 1 1 7 8Q 

l.lUjll . o u 


T a _ti nrm 
±cx iilji in 


LVG 




HV 


SN 2002cf 


—0.75,15.95 


Ia-91bg 




CL 




SN 9fin9rk 

Oli __UV.J__il_.r_ 


3.64 


Ta Tij^T"Tn 

Id 11U1 111 




CN 


N 


SN 2fin2rr 


—6.78 


T a -Ti nrm 

± cX llLJl 111 








PIN 9009^ 

Oil __rL/VJ_H_/i_i 


—7.76 


la n i~\ i' tti 

Id 11U1 111 








SN 2002ru 

Oil __ 1 1 V J _ V , 1 1 


— 5 28 1 6 62 


Ta -Tinvm 

Id 11U1 111 








SN 2002dh 

Oil __iUU--i>-LLJ 


9.21 


Ta -TinvTYi 

Id 11U1 111 






HV 


SN 2002dp 

Oil __ UVJ-jlIL 


—0 32 8 37 


Ta-nrivm 

Id 11U1 111 


HVG 


CL 


HV 


SN 9009rlf 

Oil ZrUU-i<-ll 


6.55 


la vi *~iv~m 
Id 11U1 111 








Oil __ UUiUJ 


—7.98 


Ta -Tinvm 

Id 11U1 111 








SN 9fin9rllr 
Oli _iUVJ_i<_lr_ 


— 1.23 


Ta Q1 ho- 
ld ij ± iJH 




CL 




SN 2002dn 

Oil _- V J _ V 1 1 J 


15.55 


Ta-Tirivm 

Id 11U1 111 








SN 90fl9ph 

Oil __rUU-iV_-U 


1.68 


la vi i~\ i' VTi 
Id 11U1 111 






N 


sn 2nn2ef 

Oil ^UU _i^_.l 


4.70 


Ta -Tinvm 

Id 11U1 111 




BL 


N 


SN 9n09r>h 

Oli __UVJ__iV_-ll 


6.88 


la viAyvn 
la llUl 111 








sn 2nn2ei 

Oil _- V J _ v / 1 


11.82 


Ta-Tirivm 

Id 11U1 111 








SN 2002er 


—4 58 5 26 


I £i- norm 






HV 


SN 2002et 


11.92 


Ici- norm 








SN 9fin9on 

Oil __iUl.'__L, 11 


u . uu ; y.oo 


la vh~it"wi 

Id 11U1 111 


HVG? 


CL 


N 


SN 2002fh 

Oil _- ' ' V ' — 1 I ' 


n QS 1 S fiO 
xj . c? o • i o . \j vj 


Id rj 1 t-j 




CL 




SN 9n09flr 

Oil _jUV-j1_- 


7.74 


Ta T"if~iV"m 

Id 11U1 111 






N 


SN 2002ha 

Oil _-j _-j 1 1 d 




Ta -Tinvm 

Id 11U1 111 


LVG 


BL 


N 


SN 9n09hH 

Oli __iUiJ--i-lt-- 


48 1 9 79 


la vu^wn 

Id 11U1 111 


HVG* 




N 


SN 2002he 


— 5 91 —1 03 29 3 22 


Id- norm 


HVG 


BL 


HV 


SN 2nfl2hu 

Oil __iUVJ--ll-l 


— 5.81 


1 d C/ is dd 








SN 2002hw 


—6.27 


Id- norm 








SN 2002ip- 

Oil iiUUiiJg, 


10.11 


Ta -Tinvm 

Id 11U1 111 








SN 2002iv 

Oil £j\J\J £l J V 


11.86 


Id-norm 








SN 2002kf 

Oil i- L/ \J _ I V 1 


6.81 


Ta-Tirivm 

Id 11U1 111 






N 


sn 2003n 

Oil Zj UUUJ-' 


9.98 


Ta-Tinvm 

Id 11U1 111 






N 


SN 2003K 

Oli iUlWlv 


13.43 


Ia-91T 








sn 2003TI 

Oil ZiL/vJO \J 


—2.55 


Ta-Tirivm 

Id 11U1 111 




BL 


N 


sn 2003W 

Oil ZjUVJcJ vv 


— 5 06 1 8 1 4 


T a -Ti nrm 

Id 11U1 111 








sn 2003Y 

Oil iiUUiJ J- 


— 1.74 


Ta-Q1 ho- 








SN 2003a i 

Oil iUUu CL1 


7.25 


Id I1UI 111 






N 


SN 2003ra 

Oil ^— \J \ > • 1 V A | 


—0.15 


la-norm 






HV 


SN 2003du 


17.61 


la-norm 








SN 2003fa 


-8.16 


Ia-99aa 








SN 2003gn 


-5.38 


la-norm 








SN 2003gt 


-5.07,17.61 


la-norm 








SN 2003he 


2.71,8.54 


la-norm 


LVG 


BL 


N 


SN 2003hs 


-5.49 


la-norm 








SN 2003iv 


1.76,6.58 


la-norm 


HVG? 


CL 


N 


SN 2003kf 


-7.50 


la-norm 








SN 2004E 


5.26 


la-norm 






N 


SN 2004S 


8.26 


la-norm 






N 


SN 2004as 


-4.36 


la-norm 




BL 


HV 
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Tabic Al — Continued 



SN Name 


Phase a 


SNID 


Bcnctti 


Branch 


Wang 






(Sub)Tvne b 


TvDe c 


TvDe d 


TvDe° 


SN 20041)2- 


10.34 


I fi- norm 








SN 2004bd 


10.76 


1 3,- norm 








SN 2004bk 


6.13 


± cl' lil 






HV 


SN 2004bl 


4 61 10 38 


I fi- norm 


HVG? 


CN 


N 


SN 2004br 


3.50 


± el c/ cxel 








SN 2004bv 


—7.06 9.77 


Ia-91T 








SN 2004bw 


— 10 03 6 59 


±eL 11LJ1 111 






N 


SN 2004dr 


— 6 46 1 38 18 00 


lei-norm 


HVG 


BL 


HV 


SN 2004ef 


— 5 52 8 05 


X CX 1 Awl 111 






HV 


SN 2004eo 

kJl> ZjIXVJt:i_ U 


— 5 57 1 3 1 9 


Ta.Tinnn 
xex iilji ill 








SN 2004ev 


— 7 58 1 8 80 


XeX 11LJ1 111 








SN 2004fn 

k_/ 1 N ___ LJUrrX Li. 


—2 65 2 43 


XeX 11LJ1 111 


HVG? 


BL 


HV 


SN 2004fV 


— 5 1 8 1 7 56 

J. 1U)1 1 . JU 


T a _n riT7n 
xex iilji ill 








SN 2004p-s 


0.44 


xeL iilji ill 




CL 


N 


SN 20040-u 

Oil -iiL/Urtg Li 


—4.65 


T a _ti ni'm 

X CX XXI 








SN 900^ A 


5.55 


la 








SN 200 r iM 

O i 1 U K.J i. J IV i. 


— 1 41 8 23 9 23 


Tn-Tinvtn 
xex xxljx xxi 


HVG 




N 


SN 2005W 

DIN iUUJ VV 


0.59 


T T1AVTT1 
Itl 11U1 lil 




BL 


N 


SN 2005ap- 


0.53 


Xex XXLJX XXI 






N 


SN 2005am 

ox i -_j u\j!i ei iii 


4.47 6.37 


T a _ti m'm 

XeX XXLJX XXI 


FAINT 


BL 


N 


SN 2005ao 


— 1 29 52 


la 




SS 




SN 2005hc 


1 55 7 37 


In _ t\ nviTi 
xex xxljx xxi 


LVG 


CL 


N 


SN 2005bc 

o x i __j u ^.7 ' J u I.. 


1 96 1 6 71 


Xex XXLJX XXI 








SN 2005bl 


18.07 


Ia-91bg 








SN 2005cf 


— 10.94,-2.11,-1.19,18.69 


I a,- norm 


HVG 


CN 


N 


PIN 900 

Oil ZiUUi J' 


— 7^ 10 10 


itl 11U1 ill 


HVG 


BL 


N 


SN 2005dm 

OX> ZiU\J<Jt-l.Xll 


5.23 


Ta-91 bp- 

xex j x [_/ bi 








SN 2005dv 

Oil -jUUUUV 


— 0.57 


xex XXLJX XXI 




BL 


HV 


SN 2005el 

Oil iiUVJULl 


— 6 70 1 22 8 09 


Ta.Tinnn 
xex xxljx xxi 


LVG 


CN 


N 


pin 9nn r .pr 

OlN -.uUJci 


— 9fi 1 fi7 ^ fid 


Id J X U ti 


HVG? 


CL 




Oil iUWi )H j 


— fi 01 —2 Q8 fifi 


X ex c7 exel 


HVG 






sn 9nn c .f>w 

Ol i ZiUl.'i )L, \\ 


18.23 


X ex c7 U cxet 








SN 2005cn 

Oil ^VJK.Ji. f\- 1.1 




Tn-Tinvtn 
xex xxljx xxi 








SN 2005hi 

Oil ZjIXVJ-J lij 


7.51 


Ta-Tinvm 
xex xxljx xxi 






N 


SN 2005in 


— 5.86 


XeX XXLJX XXI 








SN 900 ^Wr 


1 98 1 9 9^ 


Tn TlJ^T"Tn 
Xex ilUi ill 








SN 2005kc 


7 80 9 83 1 4 79 


xex zj x (x t-. 








SN 900 ^ki 

Oil Zj L/VJUiS-l 


1 (19 fi ^ 
A. . Ui jO .OiJ 


1 1» vi vth 

itl 11U1 ill 


LVG 


BL 


N 


SN 2005mc 

OA 1 -.UUtJlilL 


6.64 


Tn -Vinvvn 
xex xxljx xxi 






N 


SN 900^1? 


0.58 


Tn TiriVTTl 
id 11U1 ill 




BL 


N 


SN 2005ms 

Ul i iiU VJvJ lllo 


— 1 88 1 4 fi2 


Tn-Virivvn 

XeX XXLJX XXI 


HVG 




N 


SN 900^11^ 

Oli iilxUUllex 


0^1 0^ 1 7 4Q 


1 1» vi rm 
Xex 11U1 ill 


HVG 




N 


SN 2006L) 


S 70 1 (S 70 


Tn -Vinvvn 
xex xxljx xxi 




BL 


N 


SN 2006H 

Oil iUUUll 


7.01 


Tfi-Q1 ho- 

xex zj x i_/ bi 








SN 2006N 

Oil -.UUUli 


_1 sq _n 90 1 1 91 


Ta-Tinnn 

XeX XXLJX XXI 


HVG 


BL 


N 


SN 2006S 

oil -juuuu 


— 3 93 2 99 18 45 


Ia-99aa 


LVG 






SN 2006X 

Oil -jUUUiV 


3.15 


la- norm 




BL 


HV 


SN 2006ar 

Oil -jUV V ) <X\- 


7.96 


Tn - v 1 1 1 1 ■ 1 1 1 

in I1UI 111 






HV 


SN 2006ak 


8.43 


la- norm 






N 


SN 2006ax 

Oil i- \ } \ J \ 1 1 1 y\ 


— 10.07 


Ta-T1<~11"TTI 
Id I1U1 111 








SN 2006bq 


6.97,14.55,14.64 


la-norm 


HVG* 




HV 


SN 2006br 


10.62 


la-norm 








SN 2006br 

Oil iUUvL/U 


— 5 30 —4 53 2 27 


la- norm 


HVG 


CL 


N 


SN 2006bu 


4.22 


la-norm 




CN 


N 


SN 2006bw 


8.90 


la-norm 






N 


SN 2006bz 


-2.44 


Ia-91bg 




CL 




SN 2006cc 


17.67 


la-norm 








SN 2006cf 


6.28,11.09,18.69 


la-norm 






N 


SN 2006cj 


3.43 


la-norm 




SS 


N 


SN 2006cm 


-1.15,6.77 


la-norm 






N 


SN 2006cp 


-5.30 


la-norm 








SN 2006cq 


2.00 


la-norm 






N 


SN 2006cs 


2.28 


Ia-91bg 




CL 




SN 2006cz 


1.12 


Ia-99aa 
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Tabic Al — Continued 



SN Name 


Phase 3, 


SNID 


Bcnctti 


Branch 


Wang 






(Sub)Tvpe b 


Tvpe c 


Type d 


Type 


SN 2nflfidin 

k_/l^ Z.UUUU-.111 


— 7 qn s 73 1 4 fii 

1 . C»U 5 O . 1 ij , irt .U 1 


Tn-iinvtri 

Id IlUi 111 


HVG 




N 


SN 2006ef 


3.20 


Iflj-iiorm 




BL 


HV 


sn 2nnfipT 


— 8.70 


la iiui in 








SN 200fiei 


— 3 70 5 09 


I a,- norm 


HVG 


BL 


HV 


SN 2nnfiem 

k__>l> Z UUU^_ 111 


4.16 


1 <Jb rj _1 Up, 








SN 2006en 


8.55 


I a,- norm 






N 


SN 2006eu 


10.17,16.02 


I a- norm 


HVG 






SN 200fiet 


3.29 9.14 


I a- norm 


LVG 


SS 


N 


sn 2nnfipv 


1 54 1 fi 3fi 

1U. Jt:,1U.iJU 


Id 11U1 111 


HVG? 






SN 2006gt 


3.08 


Ia-91bg 








SN 2finfip-i 


4.70 


T a _ti Arm 

1(1 11U1 111 




CL 


N 


Q1N ZUUU-tY'_ 




l a Q1 ko- 
la j J- U fc-j 


HVG? 






SN 200fikf 

kJl^ ZUUU1V1 


— « 9fi — S 05 1 7 37 


Id 11U1 111 




CL 


N 


SN 200file 

1__>1>I ZjUUUIL 


— 8 fi9 1 1 Q2 

O . U C7 , J_ J. . Z 


Id 11U1 111 








SN 200filf 


— fi 30 1 4 29 


Tn -Tinvm 

Id 11U1 111 








Q1N ZUUU1I1U 


12.46 


la ilUi 111 








SN 2f)nfimn 

kJli ZjUUUIIIJJ 


5.66 


X d C/ tLv dd 








SN 2006or 


—2.79 4.93 


I a- norm 




BL 


N 


SN 2006os 


8.61 


I a- norm 






N 


SN 2006no 

k^i ^ z u uuv^u 


— 11.08 


la 








SN 200fisr 


— 2 34 2 69 

Z ■ U'i , Z . U C? 


T a _n nrm 
la 11U1 111 


HVG 


BL 


HV 


sn 2nn7A 

kJli ZjUU 1 11 


2 37 1 5 07 

Z . U 1 ,1U.U| 


Tn-iinvtri 
ia nui in 


LVG? 


CN 


N 


SN 2fin7F 


— 9 35 3 23 


Id 11U1 111 




SS 


N 


SN 9007N 

DIN ^uU I IN 


0.44 


To Q1 ho- 
ld ij 1 




CL 




SN 9007O 

kji> .zuu I v_y 


—0.33 


T i-t _n riTTn 

Id 11U1 111 




SS 


N 


SN 2007S 

k-Jl^i Zj U U i k_? 


5.18 


I a- norm 






N 


SN 2007af 

Dl'i ZUU 1 dl 


— 1 25 2 84 3 81 

_1 . Z U - Z iUtt.J.OI 


Id 11U1 111 


HVG 


BL 


N 


SN 2007a i 

Ljl>l ZUU 1 CLJ 


10.75 


T a _T"i CWTY\ 

Id 11U1 111 








SN 2007a 1 

kJl^ ZiUU 1 CL1 


3.39 


Id c/ 1 Up, 








SN 2007ar> 

k_J 1 1 iUU 1 clj_i 


9.37 


I a- norm 






N 


SN 2007au 

kJJ. 1 ZUU 1 dLl 


16.11 


Id 11U1 111 








SN 2007ax 


14.93 


T?i-Q1 hp- 

Id ij 1 [_/ bl 








SN 2007ba 


2 1 4 5 1 8 8 06 

Zj. .L^T^U. 1U,U.UU 


Ia-91bg 








SN 9007hr 

DIN ZrUU 1 


0.61 


Id IIUI III 




CL 


N 


SN 2007br1 

k_) 1 \ ZUU I l / 1 1 


7Q Q 7(1 


Id IIUI 111 






N 


SN 9n07hi 

DIN ZUU 1 Uj 


14.25 


Ta "n/^T"Tn 

Id IIUI III 








SN 9007hm 

kJl^ ZiUU i Ulll 


— 7 7Q 1 ^ OS 1 Q Qfi 


T a _n riTin 

Id 11U1 111 








SN 2007bz 


1.65 


I a- norm 






HV 


SN 2nfl7ra 

L)l "1 ZUU J t^d 


— 11 1 4 1 fi 4fi 

±1_. _lrr.iU.TtU 


Ta - ii nrm 

Id IIUI 111 








SN 9n07pcr 

DIN ZUU I 


16.17 


Ta TlJ^T"TTl 
Id IIUI III 








SN 9007ri 

kJli ZjUU 1 IjI 


_fi c;7 — 1 71 la QQ 


Ta -"nnvm 

Id 11U1 111 




CL 


N 


SN 2007cn 

kJli ZiUU 1 UU 


—4 09 85 9 55 


Ta-Tirivm 

Id IIUI 111 


LVG 


BL 


N 


SN 2007m 


— 5 82 7 84 1 5 57 

U . OZ . i . Ort ; 1U.JI 


la 








SN 2007r=; 

kJl> Z UU 1 uo 


12.15 


in IIUI III 








SN 2007fb 

UU zjuu i iu 


1 95 14 63 


la- norm 


LVG? 




N 


SN 2007fr 

t—Ji^ Zjuu i ii 


— 5 83 —1 25 


la- norm 




CL 


N 


SN 20fl7fs 

Uii ZjUU i 1j 


5.03 


la-norm 






N 


SN 2007p-c 

kJl> ZUU 1 


12.53 


Tn - v 1 1 1 1 ■ 1 1 1 

1(1 IIUI 111 








SN 2007gi 


— 7.31,-0.35,6.61 


la,- norm 


HVG? 




HV 


SN 20f)7p-k 

kjli ZjUU 1 gjlV 


— 1 72 1 9 65 

-1 . 1 — - .1 ■ ' .UU 


Ta-n<~ii"Tn 
Id IIUI 111 


HVG? 


BL 


HV 


SN 20071H 

tjl> Z UU 1 11J 


— 1 23 1 2 53 


Tn - v i . i r i vi 
in nui ill 


FAINT 


CL 


HV 


SN 2007kk 

k_) 1 \ ZUU I IS. i \ 


7.15 


Tn - 1 1 1 1 1 ■ 1 1 1 

111 IIUI 111 






HV 


SN 20071c 

k-Jl^l ZjUU 1 1*" 


_in 31 —9 40 7 43 16 39 17 37 


la- norm 


HVG 




HV 


SN 2007sl f 


-1.23 


la-norm 




BL 


N 


SN 2007on 


-3.01,-3.00,15.45 


la-norm 




CL 


N 


SN 2007qe 


-6.54,6.23,16.00 


la-norm 






HV 


SN 2007ux 


5.59 


la-norm 






N 


SN 2008C 


15.68 


la-norm 








SN 2008Q 


6.46,19.25 


la-norm 






N 


SN 2008Z 


-2.29,9.99 


Ia-99aa 








SN 2008ar 


-8.87,2.83 


la-norm 




CN 


N 


SN 2008bt 


-1.08,10.97 


Ia-91bg 




CL 




SN 2008cl 


4.24 


la-norm 




CL 


HV 


SN 2008sis 


-6.36,-4.40,-3.42,0.49,4.40,5.38,15.37 


la-norm 




BL 


N 
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Tabic Al — Continued 



SN Name 


Phase a 


SNID 


Benetti 


Branch 


Wang 






(Sub)Type b 


Type c 


Type d 


Type c 


SN 2008dx 


2.46,7.33,10.28 


Ia-91bg 


FAINT* 


CL 




SN 2008dt 


9.27 


la-norm 






HV 


SN 2008ec 


-0.24,5.70,12.51 


la-norm 


LVG 


CL 


N 


SN 2008ei 


3.29,9.13 


la-norm 


HVG* 


BL 


HV 


SN 2008s5 h 


1.26,8.96,15.76 


la 


LVG* 






SN 2008hs 


-7.94 


la-norm 









a Phases of spectra are in rest-frame days using the heliocentric redshift and photometry reference presented in table 1 
of Silverman et al. (submitted). 

b Spectral classification using the SuperNova IDentification code (SNID; [Mondin fc Tonry|2007| taken from section 5 of 
Silverman et al. (submitted). 

c Classification based on the velocity gradient of the Si II A6355 line jBenetti et al.||2005| l. "HVG" = high velocity 
gradient; "LVG" = low velocity gradient; "FAINT" = faint/underluminous. Classifications marked with a "?" are 
uncertain since light-curve shape information is unavailable. Classifications marked with a "*" use the MLCS2k2 A 
parameter I Jha et al.|2007 1 as a proxy for Ami5. 

"^Classification based on the (pseudo-)equivalent widths of the Si II A6355 and Si II A5972 lines jBranch et al.|p009| . 
"CN" = core normal; "BL" = broad line; "CL" = cool; "SS " = shallow silicon . 

"Classification based on the velocity of the Si II A6355 line I Wang et al.|2009 i. "HV" = high velocity; "N" = normal. 
f Also known as SNF20071021-000. 
sAlso known as SNF20080514-002. 
h Also known as SNF20080909-030. 



